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Introduction 


Section 2 Combinatorics 


1 Combinatorics is concerned with the arrangement, 
classification and enumeration of finite sets of objects. 


Combinatorial problems can be described under one or 
more of the following interrelated headings. 
Existence problems: 

does there exist ...? is 1t possible to ...? 
Construction problems: 

if... exists, how can we construct it? 
Enumeration problems: 

how many ... are there? (counting problem) 
-an we list them all? (listing problem) 
Optimization problems: 

it there are several ..., which is the best? 


2 Rules of counting 
Addition rule: the number of objects in a set can be 
counted by splitting the set into disjoint subsets and 
adding together the numbers of objects in each subset. 
‘Subsets are disjoint if they have no element in 
Cone ae) 
Multiplication rule: it a counting problem can be split 
into stages, each consisting of a number of options, then 
the total number of possibilities is obtained by 
multiplying the numbers of available options at each 
stave, 
Principle of inclusion—exclusion: if a set S can be split 
into two subsets A and B, not necessarily disjoint, then 
the number of elements in S is equal to 
(no. of elements in A) + (no. of elements in B) 
~ (no. of elements common to both A and B). 


The set of elements common to both A and B is called 
the intersection of A and B, and is denoted by A 1 B. 


3. Abinary word is a string of Os and Is; a single binary 
digit (O or 1) is called a bit. 


4 Analkane is a molecule with formula C,,H>,,., each 
carbon atom (C) is bonded to exactly four atoms, and 
each hydrogen atom (H) is bonded to exactly one atom 
(a carbon atom). If an alkane has four or more carbon 
atoms, then there exist different alkanes (isomers) with 
the same formula. 


5 Ann x n latin square is a square arrangement of 
n letters such that each letter appears exactly once in 
each row and exactly once in each column. 


Section 3 Representing situations 
diagrammatically 


1 A graph is a diagram consisting of points, called 
vertices, joined together by lines, called edges; each edge 
joins exactly two vertices. 

A bipartite graph is a graph in which the vertices split 


naturally into two sets, with each edge joining a vertex in - 


one set to a vertex in the other. 

In a weighted graph, each edge has a number 
associated with it called a weight (or sometimes a cost). 
A weight is a fixed number. 


2 Adigraph is a diagram consisting of points, called 
vertices, joined together by directed lines, called arcs; 
each arc joins exactly two vertices. 

In a weighted digraph, each arc has a number 
associated with it called a weight (or sometimes a cost). 


3 A network is a graph or digraph which carries some 
numerical information. This information may consist of 
weights associated with the edges or arcs, and possibly 
with some of the vertices. 


Section 4 Algorithms 


1 An algorithm is a systematic step-by-step procedure 
for solving a problem. It consists of: 

e a description of appropriate input data; 

e a finite, ordered list ot instructions; 

al SOE ER a), | 

e a description of appropriate output data. 

A heuristic algorithm is an algorithm which does not 
necessarily give a correct answer to a problem, but 
usually gives a good approximation to one. 


In a greedy algorithm, the ‘greediest’, or best, choice is 
made at each step, Bo aites of the subsequent eftect of 
that choice. 

The exhaustion algorithm examines all possible cases, 
and then selects the best. 


2 Minimum connector problem 

Given a weighted graph, find a minimum connector — a set 
of edges of minimum total weight connecting some 
specified vertices. 


Kruskal’s ereedy algorithm 7 
START with a finite set of vertices, where each pair is 
joined by a weighted edge. 

Srer 1 List the weights in ascending order. 

Step 2 Draw the vertices and weighted edge 
corresponding to the first weight in the list, 
provided that, in doing so, no cycle is formed. 
Delete the weight from the list. 


Repeat Step 2 until all | the vertices are connected, then 
eTOr 

The weighted graph obtained is a minimum connector, 
and the sum of the weights on its edges is the total weight 
of the minimum connector. 


Prim’s greedy algorithm 

Starr with a finite set of vertices, where each pair is 
joined by a weighted edge. 

Sil 
Sle 


Choose and draw any vertex. 

Find the edge of least weight joining a drawn 
vertex to one not currently drawn. Draw this 
weighted edge and the corresponding new 
verlex. 

Repeat Step 2 until all the vertices are connected, then 
STOP. 

The weighted graph obtained is a minimum connector, 
and the sum of the weights on its edges is the total weight 
of the minimum connector. 


Mh ea te ee ee ee 


3 Travelling salesman problem 

Given a weighted graph, find a route of minimum total 
weight which passes through each vertex and returns to 
the starting point. 


Heuristic algorithm for the travelling salesman 
problem 

START with a finite set of vertices, where each pair is 
joined by a weighted edge. 


STEP 1 Choose any vertex and find the edge of least 
weight emerging from it. Draw the 
corresponding two vertices and join them by two 
edges, to form a cycle. 

STEP 2 Find the edge e of least weight joining a drawn 


vertex v to a vertex w not currently drawn. Draw 
the new vertex w on the edge of the cycle 
immediately after v, moving in a clockwise 
direction. 


Repeat Step 2 until all the vertices appear in the cycle, 
then assign the appropriate weight to each edge, then 
SPOP. 

The weighted cycle obtained is a cycle through all the 
vertices, and its total weight (given by the sum of the 
weights on its edges) is at most twice the weight of the 
required minimum cycle. 


4 Efficiency of algorithms 
The efficiency of an algorithm is a measure of the time it 
takes to solve a problem. 
A polynomial-time algorithm is an algorithm for which 
the time taken is proportional to n*, for some fixed 
number k >1, or to some function f(n), where there is 
some fixed number k > 1 such that f(n) < n* for all n> N, 
for some fixed number N, where n is a parameter 
associated with the problem. Such an algorithm is 
usually efficient. Problems for which polynomial-time 
algorithms exist are said to be tractable. Problems that 
are not tractable are said to be intractable. 
An exponential-time algorithm is an algorithm for 
which the time taken is proportional to k", for some fixed 
number k>1, or to some function f(n), where there is 
some fixed number k > 1 such that f(n) > kv for all n2N, 
for some fixed number N, where n is a parameter 
associated with the problem. Such an algorithm is 
usually efficient only for problem instances of small size. 

There is an important class of problems called NP- 

complete problems. The types of problem in this class 

have the following properties: 

¢ no polynomial-time algorithm has yet been found for 
any of them; 

e no one has proved that polynomial-time algorithms 
do not exist for any of them; 

e it has been shown that, if a polynomial-time algorithm 
could be found for any one of them, then polynomial- 
time algorithms must exist for all of them; 

e it has been shown that, if it could be proved that no 
polynomial-time algorithm exists for one of them, then 
polynomial-time algorithms would not exist for any of 
them. 

It is currently believed that no polynomial-time 

algorithm exists for any of the NP-complete problems — 

that is, that NP-complete problems are intractable. 


Section 5 Mathematical modelling 


1 The mathematical modelling process 


A model is an abstraction of a real or imagined system. 

An important feature is that the important components 

of the real system should be represented by parts of the 

model. 

The mathematical modelling process may involve the 

following stages. 

Identify the practical problem. 

2 Carefully describe it in ordinary language. 

3. Formulate the description as a mathematical 
problem. 

4 Analyse and try to solve the mathematical problem. 

5 Interpret the mathematical solution in terms of the 
original problem. 

6 Examine the proposed practical solution to 
determine whether it is satisfactory. 

The mathematical formulation of a practical problem is 

called a mathematical model. 

2 Location problem 

Given a network with n vertices, representing places to be 

served by a number of facilities, the location problem 

involves one or both of the following: 

(a) finding the least number of facilities required so that 
the maximum distance between each place and its 
nearest facility does not exceed a given value d; 

(b) finding the best locations on the network for a given 
number m of facilities so that the maximum distance 
between each place and its nearest facility is a 
minimum. 


A region of a network is a set of points on the network 
such that the set of vertices of the network that can be 
reached by travelling a given distance 5 from any of the 
points in the region is always the same. A region may be 
a single point, part of an edge or parts of several edges. 


Algorithm of Christofides and Viola 

Given a location problem, START by drawing the network. 
Carry out the following steps for each distance 6. 

STEP 1 Penetrate a distance 6 from each vertex along 
the edges of the network. 

Identify the regions created, and label them 
with the combined labels of all the vertices that 
are within a distance 5 of some point in the 
region. 

Draw the bipartite graph consisting of region 
vertices on the left, representing the regions 
created in Step 2, and place vertices on the 
right, representing the vertices (places) of the 
network. Join a region vertex to a place vertex if 
the place is within a distance 6 of some point in 
the region. 

Use the bipartite graph to determine the 
minimum number and best locations of the 
facilities that can serve all the places on the 
network for the distance 58. 

In the case of location problems of type (a), carry out 
Steps 1-4 once, with 5 equal to the specified distance d. 
In the case of location problems of type (b), begin with 
5 = 0 and increase 6 in stages in such a way that each new 
region created at each stage consists of just a single 
point. STOP when the minimum number of facilities 
determined in Step 4 equals the specified number m. 


STEP 2 


STEP 3 


STEP 4 


Graphs 1 Graphs and digraphs 


Section 1 Graphs 


1 A graph consists of a set of elements called vertices, 
and a set of elements called edges. Each edge is 
associated with two vertices, and is said to join them. 
Two or more edges joining the same pair of vertices are 
multiple edges. An edge joining a vertex to itself is a 
loop. 

A graph with no leops or multiple edges is a simple 
graph 


The vertices and w are adjacent vertices if they are 
vined bv an ed > The vertices v and w are incident 


Or ee Ge =< cs 


es, oe the edge ¢ is incident with the 
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pees and j 
rwo graphs are the same if they have the same vertices 
and edges 


The complement of a simple graph G is obtained bv 
taking the vertices of G and joining two of them 
ene. ek = are Hof joined in G. 

2 Tivo graphs G and H aredsomorphic if H can be 
Le ee - “relabelling Pee Cr, that is, if there 


oS 
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is a one-one carrespondence between the vertices of G 
and these of cod that the number of edges joining 


each pair of vertices in G is equal to the number of edges 


arearsases tne c ES. pond Ng pair of vertices in ae 


3 A subgraph of a graph eee Br api h all of whose 
v ertices are vertices of G and ali at Ww hose edges are 


edg res of & 


ee 


4 Ina graph, the degree of a vertex v is the number ot 
edges incident with 7, with each loop counted twice, and 
is denoted % 82 

The degree sequence of a oe CG is the sequence 
obtained by listing the vertex degrees of G in increasing 
order, with apes as cals 


| Ticoren iE HenionaiGns lemma | 
In any graph, the sum of all the vertex degrees is | 
| estan ‘to twice the number of edges | 


A graph ist regular if its vertices all have the same degree. 


A regular alae iS f= oe or Hones of degree r, U the 
degree of each vertex is as! 


Theorem 1.2 


| ~ : ; = 
| Let G be an r-regular graph with 7 vertices; then G has 
| exactly nr/2 edges. | 
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A null graph has no edges. The null graph with 1 vertic 
denoted by N,,; it 1s regular ot degree G. 


Ny = bd ® e 
wl r 
Ns N; 
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ae toe eee consists of a single cvcle of vertices and 
edges. The cycle grapl h with # vertices is denoted by C,,; it 


1S regular of degree 2, and has 1: edges. 


in 
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icosahedron dodecahedron 


The five Platonic graphs are obtai ned by taking the 
vertices Platonic solid as the vertices 


and edges of a1 egula grapn. 


and edges oat eac oa 


tetrahedron 


dodecahedron 


icosahedron 


The Petersen graph is a 3-regular graph with 10 vertices 
and 15 edges. 


eT 


A bipartite graph is a graph whose set of vertices can be 
split into two sets A and B in such a way that each edge 
joins a vertex in A to a vertex in B. 


A complete bipartite graph is a bipartite graph in which 
each (black) vertex in A is joined to each (white) vertex in 
B by just one edge. The complete bipartite graph with r 
black vertices and s white vertices is denoted by K,.,; it 
has ¢ vertices of degree s ands vertices of degree r, and rs 


The k-cube (k-dimensional cube) Q, is the graph 
constructed by taking as vertices all 2 sequences of 0s 
and Is of length k, and joining two vertices whenever the 
corresponding sequences ditfer in just one place; it has 
2k vertices, k.2k~? edges and is regular of degree k. 


100. 101 


6 Walks, trails and paths in a graph 

A walk of length k in a graph is a succession of k edges 
of the form uv, vw, wx, ..., yz. We denote this walk by 
uvwx...yz, and refer to itas a walk between u and z. 

A closed walk in a graph is a succession of edges of the 
form uv, VW, WX, ..., YZ, ZU. | . 


A trail is a walk in which all the edges (but not necessarily 
all the vertices) are different. 

A closed trail is a closed walk in which all the edges are 
different. 

An open trail is a trail whose ends do not coincide. 

A path is a trail in which all the vertices are different. 


A cycle is a closed trail in which all the intermediate 
vertices are different. A cycle of length 3 is called a 
triangle. Any cycle of a bipartite graph has even length, 
and any graph in which every cycle has even length is 
bipartite. 


7 Agraph is connected if there is a path between each 
pair of vertices, and is disconnected otherwise. Every 
disconnected graph can be split up into a number of 
connected subgraphs, called components. 


8 A tree is a connected graph with no cycles. A tree with 
n vertices has »- 1 edges, and there is just one path 
between each pair of vertices. A tree is bipartite. 

A path graph is a tree in which there is a path through all 
its vertices. The path graph with n vertices is denoted by 


ee 
e ee e —-e—-® ee —-e—_#6 
P, P, P; P, 
Ps P, 


Section 2 Eulerian and Hamiltonian 
graphs 


1 Eulerian and semi-Eulerian graphs 

A connected graph is Eulerian if it contains a closed trail 
which includes every edge; such a trail is an Eulerian 
trail. 

A connected graph is semi-Eulerian if it contains an 
open trail which includes, every edge; such a trail is a 
semi-Eulerian trail. 


Theorem 2.1 | 
| A connected graph is Eulerian if and only if each | 
| vertex has even degree. 
ae et 


ae S Soe een eran rer eee or 


Theorem 2.2 | 
A connected graph is semi-Eulerian if and only if it — 
_ has exactly two vertices of odd degree. | 


Theorem 2.3 | 
An Eulerian graph can be split into cycles, no two of | 
which have an edge in common. 


| 
i} 


2 Hamiltonian and semi-Hamiltonian graphs 

A connected graph is*Hamiltonian if it contains a cycle 
that includes every vertex; such a cycle is a Hamiltonian 
cycle. 

A connected graph is seri-Hamiltonian if there is a 
path, but not a cycle, which includes every vertex; such a 
path is a semi-Hamiltonian path. 


| Theorem 2.4: Ore’s theorem 
| Let G be a simple connected graph with n vertices, 
_ where n 2 3. If deg v + deg w >n, for each pair of non- 


| adjacent vertices v and w, then G is Hamiltonian. 
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Section 3 Digraphs 


1 A digraph consists of a set of elements called vertices, 
and a set of elements called arcs. Each arc joins two 
vertices in a specified direction. 

Two or more arcs joining the same pair of vertices in the 
same direction are multiple arcs. An arc joining a vertex 
to itself is a loop. 

A digraph with no loops or multiple arcs is a simple 
digraph. 

The vertices v and w are adjacent vertices if they are 
joined (in either direction) by an arc e. An arc ¢ that is 
directed from v to w is incident from 7 and incident to 7; 
< is incident to e, and w is incident from e. 


The underlying graph of a digraph D is the graph 
obtained by replacing each arc of D by the 
corresponding undirected edge. 


Two digraphs are e the same if they have the same vertices 
and arcs. : 


2 Two digraphs C and D are isomorphic if D can be 
obtained by relabelling the vertices of C, that is, if there is 
a one-one correspondence between the vertices of C and 
those of D such that the arcs joining each pair of vertices 
in C agree in both number and direction with the arcs 


joining the corresponding pair of vertices in D. 


3 A subdigraph of a digraph D is a digraph all of whose 
vertices are vertices of D and all of whose arcs are arcs ot 
oe 


4 In a digraph, the out- Hearne of a vertex v is the 
number of arcs incident from v, and is denoted by 
outdeg v; and the in-degreeof » is the oilags of ares 
incident to v, and is denoted by indeg v. Each loop 
contributes | to. both the in- degree and Te out: degree of 
the corresponding vertex.” 

The out-degree sequence of a digraph D is the sequence 
obtained by listing the out-degrees of D in increasing 
order, with repeats as necessary. The in- degree 
sequence of D is detined analogo ously. 


5 
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Theorem 3.1: Hencenoene dilemma 


In any digraph, the sum of all the out-degrees and the | 4 
sum of all the in-degrees are each equal to the © 


number ot arcs. 
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6 Walks, trails and paths in a digraph 

A walk of length k in a digraph is a succession of k arcs 
of the ais uv, VW, WX, ..., yz. We denote this walk by 
unwx... yz, and refer to itas a walk from u to z 

A eee wail in a digraph is a succession of arcs of the 
forma 00 eae a. 

A trail is a walk in which all the arcs (but not necessarily 
all the vertices) are different. 

A closed trail is a closed walk in which all the arcs are 
different. 

A path is a trail in which all the vertices are different. 


A cycle is a closed trail in which all the intermediate 
vertices are different. 


7 A digraph is connected if its underlying graph is a 
connected graph, and is disconnected otherwise. It is 
strongly connected if there is a path between e each pair 
of vertices. 


8 Eulerian digraphs 
A connected digraph is Eulerian if it contains a closed 


trail which includes every arc; such a trail is an Eulerian 


trail. 


| ieee ep 
_ A connected digraph is Eulerian if and only if, for sacs 
| vertex, , the o out- it-degree equals the in- degree. 


- Theorem 3.3 
_ An Eulerian digraph can be split into cycles, no two of 
which have z an arc in common. 
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9 Hamiltonian digraphs 

A connected digraph is Hamiltonian if it contains a 
cycle which includes every vertex; such a cycle is a- 
Hamiltonian cycle. 


Section 4 Matrix representations 


1 Adjacency matrix of a graph 

Let G be a graph with n vertices labelled 1, 2 . The 
adjacency matrix A(G) is the 1 x n matrix in oli the 
entry in row? and column 7 is the number of edges 
joining the vertices i and ;. The matrix is symmetrical 
about the main diagonal (top left to bottom right). If the 
graph has no loops, then each entry on the main 
diagonal is 0, and the sum of the entries in any row or 
column is the degree of the vertex corresponding to that 
row or column. 


2 Adjacency matrix of a digraph 

Let D be a digraph with n vertices labelled 1, 2 ..., n. The 
adjacency matrix A(D) is the n x m matrix in ee the 
entry in row i and column j is the number of arcs from 
vertex 7 tO vertex |. The matrix is not necessarily 
symmetrical about the main diagonal. If the digraph has 
no loops, then each entry on the main diagonal is 0, and 
the sum of the entries in any row is the out: degree of the 
vertex corresponding to that row, and the sum of the 
entries in any column is the in-degree of the vertex 
corresponding to that column. 


3 Incidence matrix of a graph 
Let G be a graph without loops, with n vertices and 
m edges. The incidence matrix B(G) is the 1 x m matrix 
in which the entry in row i and column 7 is 

1 if the vertex 7 is incident with the edge /, 

Q otherwise. 


4 Incidence matrix of a digraph 
Let D be a digraph without loops, with m vertices and 
m arcs. The incidence matrix B(D) is the 1 x mm matrix in 
which the entry in row 7 and column J is 

1 if arc j is incident from vertex 1, 
-l1 ifarcyis batefersi ie to.vertex L, 

0 otherwise. 


Networks 1 Network flows 


Section 1 Connectivity 


1 Edge connectivity 
The edge connectivity A(G) of a connected graph G is the 
smallest number of edges whose removal disconnects G. 


ayy 


A cutset of a connected graph G is a set S of edges with 

the properties: 

(a) removal of all the edges in S disconnects G; 

(b) removal of some (but not all) of the edges in S does 
not disconnect G. 

Two cutsets Ble a ee need not necessarily have the 


The aes meoany x(G) of a graph G is the size of the 
smallest cutset of G. 


A single edge whose removal disconnects a 
called a bridge. 


graph 1s 


2 Vertex connectivity 

The connectivity (or vertex connectivity) K(G) of a 

connected graph G (other than a complete graph) is the 

smallest number of vertices whose removal disconnects 

G. 

The connectivity «(K,,) of the complete graph K,, (1: 

H— I. 

A vertex cutset of a connected graph G is a set S of 

vertices (not the whole set of vertices) with the foll owing 

properties: 

(a) removal of all the vertices in S disconnects G; 

(b) removal of some (but not all) of the vertices in S does 
not disconnect G. 

Two vertex cutsets of a graph need not necessarily have 

the same number of vertices. 

The connectivity x(G) of a graph G is the size of the 

smallest vertex cutset of G. 


ao) tS 


A single vertex whose removal disconnects a graph is 
called a cut vertex. 


> 


_ Theorem 1.1 : 
| For any connected graph G, : 
| K(G) < A(G) < &(G), 


! _ where d(G) is the smallest vertex « degree i in G. 


4 st-paths in a graph 
Let G be a connected graph, and let s and ft be vertices of 
G. A path between, sand ¢ is called an st-path. Two or 
more st-paths are edge- disjoint if they have no edges in 
common, and vertex- -disjoint if they have no vertices in 
common (apart from s and f). Certain edges separate s 
from t if the removal of these edges destroys all paths 
between s and ¢. Similarly, certain ee (excluding s 
and t) separate s from tif the removal of these vertices 
destroys all eee between s and t. 

| Theorem clea. Meroe & aheorem (edee fon 

Let G be a connected graph, and let s and f be vertices | 
of G. Then the maximum number of edge-disjoint sf | 

paths is equal to the minimum number of edges 

separating » s from t. i 


If we can find k edge-disjoint st-paths and k edges 
separating s from t (for the same value of k), then k is the 
maximum number of edge-disjoint st-paths and the 
minimum number ot edges separating s from f. These k 
edges separating s from ft necessarily form a cutset. It 
follows that, when looking for them, we need consider 
only cutsets whose removal disconnects G into two 
components, one containing s and the other containing 1. 
Corollary of Menger’s theorem (edge form) 

A connected graph G has edge connectivity / if and 

only if every pair of vertices in G is joined by 1 or more 

edge- disjoint paths, and at least one pair of vertices is 

jomned by exactly | edge-disjoint paths. 


5 st-paths in a digraph 
Let D be a connected digraph, and let s and t be vertices 
of D. A path trom s tot is called an st-path. Two or more 
st-paths are arc-disjointif they have no arcs in common, 
and vertex-disjoint if they have no vertices in common 
(apart from s and f). Certain arcs separate s from t if the 
removal of these arcs destroys all paths from s to ¢. 
Similarly, certain vertices (e xcluding s and t) separate s 
from ft if the removal of these vertices destroys all paths 
from s tot. 

Theorem 1.3: Menger’s theorem ae on 

Let D be a connected digraph, and let s and t be | 

vertices of D. Then the maximum number of arc- | 
_ disjoint st-paths is equal to the minimum number of | 
Parcs SI PEM IN s from ft. 


If we can find k arc-disjoint ee and k arcs separating 
s from t (for the same value of k), then k is the maximum 
number of arc-disjoint sf-paths and the minimum number 
of arcs separating s trom i 
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| THEOTein 1.4: ONeneers Benen for graphs (vertex 


form) 

Let G be a connected graph, and let s and t be non- 

adjacent vertices of G. Then the maximum number of | 
_ vertex-disjoint st-paths is equal to the minimum 
| mone of vertices Ss Sepeue bys: from tf. 


If we can find k vertex-disjoint st- ae and k vertices 
separating s s trom i (tor the same value of ky, then k is the 
maximum number of vertex-disjoint st-paths and the 
minimum number of vertices separating s from t. These k 
vertices separating s from f¢ necessarily form a vertex 
cutset. It follows that, when looking for them, we need 
consider only vertex cutsets whose removal disconnects 
G into two or more components, one containing s and the 
other containing t. 


Corollary of Menger’s theorem for graphs (vertex | 
form) 

A connected graph G (other than a complete graph) 
has vertex connectivity k if and only if every non- 
adjacent pair of vertices in G is joined by k or more 
vertex-disjoint paths, and at least one non-adjacent 
pair of vertices 1s joined by exactly k vertex-disjoint 


paths. 


en 


i | 


Theorem 1.5: Menger’s theorem for digraphs (vertex 
form) 

Let D be a connected digraph, and let s and t be non- 
_ adjacent vertices of D. Then the maximum number of | 
| vertex- disjoint st-paths is equal to the minimum | 
pee! of vertices SENSE $ from t 


8 Let G be a graph with n vertices and m edges. Then G 
has optimal connectivity if K(G) = 2m/n. 

All graphs with optimal connectivity are regular graphs, 
but not every regular graph has optimal connectivity. 


Section 2 Flows in basic networks 


1 A basic network is a connected digraph N satistving 

the conditions: 

(a) N has exactly one source and one sink; 

(b) each arc e of N is assigned a positive number c(e), 
called the capacity of e, which represents the 
maximum amount of the commodity that can flow 
through the arc in a given time. 

(Capacities are shown in bold type in diagrams.) 


2 A flow ina basic network N with source S and sink T is 
an assignment, to each arc e of N, of a non-negative 
number f(e), called the flow along the arc ¢, satisfy ‘ing: 
(a) the feasibility condition: the flow along each arc 
does not exceed the capacity of the arc, that is, 
fle) < c(e), for each arc e of N; 
(b) the flow conservation law: for each vertex V other 
than S and T, the sum of the flows along the arcs into 
V is equal to the sum of the flows along the arcs out 
of V. 
Such a flow is sometimes called a flow from S to T or an 
ay T)-flow : 


An arc ¢ is saturated if f(e) = c(e), and it is unsaturated if 


fle) < cle). 


(A saturated arc is shown by a thick line in diagrams.) 


The total flow out of the source S is always equal to the 
total flow into the sink T. This common value is called the 
value of the flow, and represents the amount of the 
commodity flowing through the network. 

A maximum flow is a flow of largest possible value. 
(There may be several different flows with the same 
largest value.) 


3 A flow-augmenting path in a basic network with 

source S and sink T is a path from S to T consisting of: 

¢ forward arcs: unsaturated arcs directed aioue the 
path; 
and, possibly, 

* backward arcs: arcs directed against the direction of 
the path and carrying a non-zero flow. 


4 ‘Transforming to a basic network 

An undirected network is one in which flow is allowed in 

either direction along each edge. A network containing 

both directed and undirected edges is a mixed network. 

The transformations to convert a network into a basic 

network are summarized below, and should be carried 

out in the order given. bes 

e Replace each undirected edge by two arcs, one in 
each direction, both with the same capacity as the 
original undirected edge. 


k 
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e—- as oD - " ape ee 
I 


e Replace each vertex V with a capacity restriction k by 
two vertices V,; and V, joined by an arc trom V, to V; of 
capacity k. All arcs directed towards V are Meera 
towards V,, and all arcs directed away trom V are 
directed away from V3. 


o Pe os 2 oe 
a a UL” a aS tA 
i ge : a k A 
Sec Paes 
ce ace ue Vi Vea 
Ed a ” nN 
“e o ‘e 


e If there are several! cone Soy be eee Om tem 
to a new super-source S. If there are sev eral sinks tee 
der roe mere chro devsras - new super-sink T. To oer 

new arc S assign a capacity equal to the sum of the 
capacities out of S; to each new arc TT assign a 
capacity equal to the sum of the capacities into T;. 
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General procedure 

A variety of network flow problems are solved by the 
following procedure. 

Step 1 Transform the network into a basic network 
using the above transformations where 
necessary. 

Solve this basic network problem. 

Interpret the solution in terms of the original 
network. 


ho 


STEP 
STEP 


wo 


Section 3 Maximum flows and minimum 
cuts 


1 Acut ina basic network with source S and sink T isa 
set of arcs whose removal separates the network into two 
disjoint parts X (containing S) and Y {containing 1). 

The capacity of a cut is the sum of the capacities of those 
arcs in the cut which are directed from X (the part 
containing S) to Y (the part containing 1). 

A minimum cut is a cut of smallest possible capacity. 


2 Flows and cuts 
An important connection between flows and cuts is: 
the value of any flow < the capacity of any cut. 


3 Maximum flows and minimum cuts 
The connection between maximum flows and minimum 
cuts is: 
the value of any maximum flow < the capacity of any 
nunimum cut. 
If we can find a flow with value k, and a cut with capacity k 
(for the same value of k), then: 
« the flow is a maximum flow; 
e the cut is a minimum cut. 
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Theorem 3.1: max-flow min-cut theorem | 

_ In any basic network, the value of a maximum flow is | 
_ equal to the capacity of a minimum cut. 


[NEMO face Mei aire He eae 


An aaporennt consequence ae ne Ghesian is aes if you 
can find a minimum cut in a network, then you can 
immediately deduce the value of a maximum flow. 


4 Maximum flow algorith 
START with any flow (possibly the zero flow). 


Part A: labelling procedure 

Start from the source S and proceed a 
allow further flow to unlabelled ver ttices, ch hosen 
alphabetically. Use a forward arc where possible, but if 
none exists try using an arc carrying a non-zero flow as a 
backward arc. 

Label each vertex with the letter denoting the previous 
vertex and a number giving the size of the change of flow 
possible from that verte 

If a vertex is reached from which no further prog 
possible, cross through the current label of that» 
and try backtracking to the vertex previously ae 
and continue from there. | 

If the sink T becomes labelled, we say we have 
breakthrough; go to Part B. 


If there is no breakthrough, stop: 
e the existing flow is a maximum flow; 


* a corresponding minimum cut separates S and the 
currently labelled vertices from the other vertices. 


Part B: flow-augmenting procedu 

Augment the flow along the path of the labelled vertices 
S,...T by the largest flow consistent with the labels. 
Remove the labels and return to Part A. 


Section 4 Networks with lower and 
upper capacities 


1 In some networks, it is necessary to specify both the 
maximum and the minimum allowable flow along each 
arc. The minimum allowable flow along an arc is the 
lower capacity of the arc, and the maximum allowable 
flow is the upper capacity of the arc. A basic network can 
be regarded as a network of this type in which all lower 
capacities are zero. 


2 A flow in a network with lower and upper capacities is 
an assignment of non-negative numbers to the arcs in 
such a way that the following conditions are satisfied: 

(a) the feasibility condition: the flow along each arc is 
not less than the lower capacity and is not more than 
the upper capacity of that arc: 

lower capacity < flow < upper capacity; 

(b) the flow conservation law: the flow into each vertex 
(other than S or T) is equal to the flow out of it. 

(We represent both flows and capacities on the same 

diagram: the flow appears between the lower capacity 

and the upper capacity, and the capacities are shown in 
bold.) 

Terms such as the value of a flow and a maximum flow 

are defined as before. 


3 If we can find a flow in a network with lower and upper 
capacities, then we can increase this flow to a maximum 
flow using the maximum flow algorithm to increase the 
flow along flow-augmenting paths, provided that when 
using a backward arc we do not decrease the flow in that arc to 
a value less than the lower capacity of the arc. 

Finding an initial flow is not always easy, since we can no 
longer choose the zero flow as the initial flow. There may 
not even be such a flow. 


4 Feasibility 

A network N with lower and upper capacities is feasible if 
there exists a flow satisfying all the capacity restrictions. 
If no such flow exists, then the network is infeasible. 

The idea behind the feasibility algorithm is to construct 
from the given network N a basic network N* in such a 
way that if there is a flow in N* which satisfies certain 
conditions (see feasibility theorem below), then the 
network N is feasible. 
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Feasibility algorithm 
Construct N* from N by the following procedure. 
STAGE 1 Take the vertices of N* to be the vertices of N, 


together with two new vertices: a new source S* 
and a new sink T*. 


S* 
® 


a 


STAGE2 Replace the lower and upper capacities | and u 
on each arc VWof N by a single capacity u —1 on 
the corresponding arc of N*. 


V Lu 


STAGE3 For each arc VW with lower capacity 1 in N, add 
in N* an arc S*W with capacity | and an arc VI* 
with capacity 1. 


S* 
e---J 
«= : 
u-l SW 
oe Tx 
oe 


STAGE4 In N*, add anew arc from T to S with infinite 
capacity. 


, N 
i \ 
- \ 


e Toshow that N is feasible, find a flow in the basic 
network N* in which all the arcs out of S* and all the 
arcs into T* are saturated. (Do this by inspection or 
by using the maximum flow algorithm.) 

Convert this maximum flow in N* to a flow in N by 
‘undoing’ the above construction. 

If required, increase the flow in N to a maximum 
flow. (See item 3.) 

e Toshow that N is infeasible, show that in N* there is 
no flow from S* to T* in which all the arcs out of S* 
and all the arcs into T* are saturated. 


Theorem 4.1: feasibility theorem 
A network N with lower and upper capacities is | 
feasible if and only if, in the related basic network N*, | 
there is a flow from S* to T* such that all the arcs out 

of S* and all the arcs into T* are saturated. ; 


5 Let N be a network with lower and upper capacities, 
and let C be a cut which separates N into two disjoint 
parts: X (containing S) and Y (containing T). Then the 
capacity of the cut Cis 

(the sum of the upper capacities of the arcs of C 

directed from X to Y) 

— (the sum of the lower capacities of the arcs of C 

directed from Y to X). 

(If the lower capacity of each arc is zero, this definition 
coincides with the definition in Section 3, item 1.) 
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Theorem 4.2: generalized max-flow min-cut 
theorem 


In any feasible network with lower and upper 
capacities, the value of a maximum flow is equal to 
the capacity of a minimum cut. 
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eiric elements 


1 The dime 


nsion of a space is the minimum number of 
coordinates needed to describe the position of a typical 
point in that space. 

Basic geometric elements in 0, 1, 2 and 3 dimensions are, 
respectively, points, line segments, plane segments and space 
segments, The simplest type of plane segment is a 
triangle and the simplest type of space segment is a 
tetrahedron: 


? cep pes 

A set S of points in a Euclidean space is convex if, for 
each possible choice of two points P and Q in S, ail the 
points on the line segment joining P and Q also hie in S. 

To show that a shape Sis not convex, it is sufficient to find 
one pair of points P,Q in S such that the line segment PQ 
does not hie wholly in S. 


3. Convex hulls 
Let S be a set of points in a Euclidean space. The convex 


hull of S is the smallest convex set containing S; it is 
denoted by op 


on 
2 
oe 
ae 
a 
cas) 
oe 
© 
et 
pp) 
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ag 
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al] lie ona line. Then the convex hull (S) is a convex 
polygon. 


| 
bea Anite set of points in the plane that do not | 


reais an object in space bounded by plane 
it consists of all the points on, and 
lane segments. 


elie ( ts faces; 

| , the boundir ‘BP pla 
Theerem 1.2 | 
Let S be a finite set of points in space that do not all he | 


in the same plane. Then the convex hull (S) is a convex | 
| polyhedron. | 


4 The points and line segments bounding a convex hull 
can be interpreted as the vertices and edges of a graph, 
called the graph of the convex hull. 


A single pow. is the onl ¥ y type of convex zero-dimensional 


set, and i is called a @«: ley. 


A oe eginent is the simplest a Ee of convex one- 
dimer anal set, and is called a 7 plex. 
A triangle is the pale est type - convex two-dimensional 


ei, mee is called a ?-simplex. 
A. tetrahedron is ne simplest type of convex three- 


dimensional set, and is called a 3-sime 


Ann-simplex, the simplest type of convex 1-dimensional 
set, is the convex hull of a set of m + 1 points, placed in the 
most general possible positions in n-dimensional space; 
each line ale between pairs of these points bounds 
x. Its o es the ee a graph Kj. ;. 


the n-simple 


ard n-arthotape ©, is constructed as follows: 
O, is the unit interva! a eae atx =QOand x =1; 
O, is the unit square with vertices at the pon 
(0, 0), (1, 0), (0, 1), Ds 
O, is the uiit cube with vertices at the points 
(0;.0;,0),.00) 07 1))(0; 1.0), (07.45. 2), 
(1,0;05,.(1, 054 ),0,4,0), Cl, 1): 


In general, O,, has 2" vertices, each with n coordinates, 
each equal to 0 or 1. The graph of O,, is the n-cube. 


7 The n-dimensional analogue of a polygon or 
polyhedron is a peiytope. 


z 


A cross pelytope is a convex polytope formed by taking 
the convex hull of the ‘cross’ formed by going a certain 
distance in a positive and a negative direction along each 
coordinate axis. The standard cross polytope is 
constructed by going out a distance 1 along each axis. 

In n-dimensional space, a standard n-dimensional cross 
polytope is obtained. 

The graph of a cross polytope is formed by drawing an 
edge between each pair of points except those on the 
same coordinate axis. In the one-dimensional case there 
are just two vertices, corresponding to the two points, and 
no edges; so the graph is the null graph N>. In the two- 
dimensional case a square with four vertices each of 
degree 2 is obtained; so the graph is the cycle graph C, 
In the three-dimensional case an octahedron with six 
vertices each of degree 4 is obtained, which can be drawn 
as the standard octahedron graph. The graph of an 
n-dimensional cross polytope is sometimes called an 
n-octahedron. 


away from a given point: 
a 1-sphere is a line segment; 
a 2-sphere is a disc; 
a 3-sphere is a solid sphere (in the usual sense of the 
word). 


1 A packing of shapes or objects is constructed from 
several copies of a single shape or object arranged to 
give a regular repeating two-dimensional pattern. If a dot 
is placed at the centre of each shape or object, a regular 
array of points, called a point lattice, is obtained. 
Packings of shapes in the plane include the triangular, 
square and hexagonal packing arrangements. 


2 Atil 
gaps in it 
A tiling (or tess 
whole plane with t tile s in 
gaps or overlaps. 


> is a planar shape of Hite area with no holes or 


ellation) o of the plane is a covering of the 
in such a way that there are no 


Two planar shapes are congruent if they are the same 
size and shape, that is, if one can be superimposed 
exactly on the other, possibly after turning one over. A 
given shape tiles the plane if the plane can be tiled using 
only congruent copies of that shape. 


A tiling with polygons 1s an e¢ ee tent edge tiling if each 
edge of each polygon is an edge of another aes and 
no two polygons meet along more than one edge. 


egular polygen is a polygon in which all the edges are 
of equal length, a the interior angles between each 
pair of adjacent edges are all equal. A regular polygon is 
necessarily convex, but a convex polygon need not be 
regular. 

The following table gives the interior angle @ in radians 
and degrees of an n-sided regular polygon. 


4 5 6 7% 8 9 0 12 


Q t/3- n/2- 3x /5 Qn /3 Sit /7 3n /4 In /9 An /5 51 /6 
60° 90° 108° 120° 129° 135° 140° 144° 150° 
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In any edge-to-edge tiling using regular polygons, the 
tiles fit together in such a way that each edge belongs to 
just two tiles and each vertex belongs to at least three 
tiles. 

A vertex type is the arrangement of regular polygons 
around a vertex 


3.3.3.3.6 


3 Regular tilings 
A regular tiling is an edge-to-edge tiling with congruent 
regular polygons. 
There are only three tvpes of regular tiling; with triangles, 
squares and hexagons: 

triangular, vertex type: 3.3.3.3.3.3 

square, vertex tvpe: 4.4.4.4 

hexagonal, vertex type: 6.6.6 


3.3.3.4.4 


4 If two or more types of regular polygon are present at 
each vertex, there are a further eighteen possible vertex 
types. 

A semi-regular tiling is an edge-to-edge tiling with 
regular polygons, not all congruent, such that each vertex 
tvpe is the same (up to reflection). There are eight semi- 
regular tilings: 


(c) 3.3.4.3.4 (d) 3.4.6.4 


(f) 3.12.12 
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5 A demi-regular tiling is an edge-to-edge tiling with 
regular polygons in which more than one vertex type 
occurs. There are infinitely many demi-regular tilings. 


6 Dual tilings 

Edge-adjacent polygons are polygons that share a 
common edge. 

Given an edge-to-edge tiling with regular polygons, a 
dual tiling is constructed by joining the centre of each 
polygon to the centres of its edge-adjacent polygons, 
using line segments. 


Tne Laves tilings are the duals of the regular and semi- 
regular tilings. Three of them are regular tilings — the 
triangular and hexagonal tilings are duals of each other, 
and the square tiling is self-dual. 

The duals of the eight semi-regular tilings are not 
regular, semi-regular or demi-regular. 


7 A polygonal animal is a connected arrangemeni of 
non-overlapping regular polygons in edge-to-edge 
contact. If it comprises exactly 1 polygons, it is called an 
n-animal. 


formed from — 


equilateral triangles n-iamond 
squares n-Omino 
regular hexagons n-hex 


The graph of a polygonal animal is formed by placing 
vertices in the centres of the polygons and joining the 
vertices in edge-adjacent polygons. 


Section 3 Tilings at the Alhambra 


1 Arotation of order k is a rotation through 21/k radians 


or 360/ k degrees. 
CC ¥. oe P . 
Section 4 Polyhedra 


1 Regular polyhedra — 

A regular polyhedron is a convex polyhedron in which ail 
the polygonal faces are congruent regular polygons, and 
in which each vertex has exactly the same arrangement 
of polygons around it. 


There are exactly five regular polyhedra, known as the 
Piatonic solids. (See also Graphs 1, Section 1, item 5.) 


polyhedron vertices ed ges faces 
tetrahedron | 4 6 4 
octahedron | 6 12 8 
cube 8 12 6 
icosahedron 12 30 20 


dodecahedron 20 30 12 


The graphs of the regular polvhedra are called the 
Platonic graphs. | 


2 Dual polyhedra 

The dual of a convex polyhedron is constructed by 
placing a new vertex at the centre of each face of the 
original polyhedron, and joining a pair of new vertices 
with a line segment whenever the corresponding faces of 
the original polyhedron are edge-adjacent; thus the roles 
of the vertices and faces are exchanged. 

The duals of all the Platonic solids are Platonic solids. 


3 Semi-regular polyhedra 

A semi-regular polyhedron is a convex polyhedron in 
which all the polygonal faces are regular polygons, not all 
congruent, such that each vertex has the same 
arrangement of polys 
There are intinitel 
comprising two infinite sets (the prisms and antiprisits ) 
and one finite set of thirteen polyhedra called the 
Archimedean solids. 


A prism is constructed by taking two parallel congruent 
regular polygons, aligned one above the other, and 
joining, corresponding vertices with new edges of length 
equal to that of the sides of the polygons. This results in 
an ‘equatorial strip’ of squares with a polvgonal top and 
base. 

An antiprism is constructed similarly from two parallel 
congruent regular polygons. The polygons are not 
aligned; the vertices of the upper polygon are placed 
above the midpoints of the edges of the lower one. The 
vertices are Joined to produce an ‘equatorial strip’ of 
equilateral triangles. 


4 Results about pelyhedra 


' Theorem 4.1: Euler’s polyhedron formula | 
Let v,e and f denote, respectively, the numbers of | 
vertices, edges and faces of a convex polyhedron. 

Then 


| 5 | 


The degree of a face of a polyhedron is the number of 
edges around it. 


Theorem 4.2: handshaking lemma for polyhedra 
In any polyhedron, the sum of all the-face degrees is — 
qual to twice the number of edges. : 


Fans 


Theorem 4.3 

There are only five regular polyhedra: 

e three with triangular faces — the tetrahedron, the 
octahedron and the icosahedron; 

* one with square faces — the cube; 


' 
| 
| 
| 
f 
| 


— ® one with pentagonal faces — the dodecahedron. 


Section 5 Incidence structures 


1 An incidence structure 5(P, L) consists of two sets P 
and L, of different types of object, and an incidence 
relation that shows, for each pair of objects, one from P 
and one from L, whether or not these two objects are 
incident. 

We often describe the elements of the two sets P and L 
as ‘points’ and ‘lines’, respectively. 


The degree of a ‘point’ 1s the number of ‘lines’ incident 
with it, and the degree of a ‘line’ is the number of ‘points’ 
incident with it. 


An incidence structure is regular if every ‘point’ has the 


be 


same degree d, and every ‘line’ has the same degree dj). 


2 nepresenting incidence structures 

Let S(P, L) be an incidence structure, with 1,, ‘points’ 
labelled 1 to m, and mn; ‘lines’ labelled | to 1). The 
incidence matrix B(S) of Sis the 1, x 7) matrix in which 
the entry in row 7 and column j is | if the ‘point’ labelled 7 
is incident with the ‘line’ labelled 7, and 0 otherwise. 


Let A be any matrix with m7 rows and nm columns. Then the 
transpose of A, denoted by A‘, is the matrix with 1 rows 
and m columns obtained by interchanging the rows and 
columns of A; that is, the entry in row 1 and column j of A 
is the same as the entry in row] and column / ot A’ 

3 Dual incidence structures 

Let S(P, L) be an incidence structure. Then the incidence 
structure S*(L, P) obtained by exchanging the roles of P 
and L, is called the dual of S(P, L). 


Theorem 5.1 

Let S be any incidence structure, and S* its dual. Then: 

(a) the incidence matrix of S* is the transpose of that 
or 

(b) the sum of the degrees of the ‘points’ of 5 is equal 

to the sum of the degrees of its ‘lines’. 


4 An incidence structure can be represented by a block 
table consisting of a set of columns, each of which 
represents a ‘line’ by a list of those ‘points’ with which it is 
incident. 


5 Finite projective geometry 


A finite projective plane is a regular incidence structure 

for which: 

a) for each pair of distinct 'points', there is exactly one 
‘line’ incident with both; 

b) for each pair of ‘lines’ there is exactly one ‘point’ 
incident with both. 

The Fano plane is a finite projective plane with seven 

points and seven lines, each line being incident with three 


points, and each point with three lines. 


re 


Graphs 2 Trees 


Section 1 Tree structures 


1 A tree is a connected ue that has no cycles. 


| Theorem 1.1: Seavaienl ahaa of a tree 

| Let T be a graph with n vertices. Then the following 
_ statements are equivalent. | 
| ¢ Tis connected and has no cycles. | 
¢  Thas n—-1 edges and has no cycles. ! 


| ° “is connected and has n — 1 edges. 
| a 
, ® Tis connected and the removal of any edge 


disconnects T’. 


oh 
IS 


| e Any two vertices of T are connected by exactly 
one path. 
e fcr 


ains no cycles, but the addition of any new 


creates a cycle. 


A rooted tree (or branching tree) is the hierarchical 
structure in which one vertex is singled out as the starting 
point, and the branches fan out from this vertex. The 
starting vertex is the root, and is often represented by a 
small square. 

A sorting tree is a branching tree that arises when a 
succession of choices is made, each dependent on the 
previous one. 


Section 2 Counting trees 


1 Counting labelled trees 

There is a one-one correspondence between labelled 
trees with n vertices and Priifer sequences — sequences 
of 1 - 2 numbers (aj, ao, -.., a2), where each a; is one of 
the integers 1, ..., 1 (allowing repetition). 


To construct a Prifer sequence from a given labelled 
tree 
Step 1 Find the vertices of degree 1 and choose the one 


with the smallest label. 

Look at the vertex adjacent to the one just 
chosen, and place its label in the first available 
position in the Priifer sequence. 


STEP 2 


SteP3 Remove the vertex chosen in Step 1 and its 
incident edge, leaving a smaller tree. 

Repeat Steps 1-3 for the remaining tree, continuing until 

there are only two vertices left, then STOP: the required 

Priifer sequence has been constructed. 


To construct a labelled tree from a given Prifer 

sequence containing n — 2 numbers 

STEP 1 Draw n vertices, labelling them from 1 to 1, and 
make a list of the integers from | to n. 

Step 2 Find the smallest number which is in the list but 
not in the Priifer sequence, and also tind the first 
number in the sequence; then add an edge 
joining the vertices with these labels. 

SreP 3. Remove the first number found in Step 2 from the 
list and the second number found in Step 2 trom 
the sequence, leaving a smaller list and a smaller 
sequence. | 

Repeat Steps 2 and 3 for the remaining lst and 

sequence, continuing until there are only two terms left 

in the list. Then join the vertices with these labels and 

STOP: the required labelled tree has been constructed. 


Theorem 2.1: Cayley’s theorem 


The number of f labelled b trees with 1 vertices is 1” ~ | 
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2 ‘Counting binary trees 

A binary tree is a rooted tree in which the number of 
edges descending from each vertex is at most 2, and a 
distinction is made between left-hand and right-hand 
branches. 

Let u,, denote the number of binary trees with n vertices. 
Then u, = 1, u.=2, and for n 23, 


cri em YY : ; fag A whew ¢ y 
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3 Central and bicentral trees 

To find the ‘middle’ of a tree, remove all the vertices of 
degree | together with their incident edges, and repeat 
this process until either a single vertex, called the centre, 
or two vertices joined by an edge, called the bicentre, is 
obtained. A tree with a centre is called a central tree, and 
a tree with a a iS eee a bicentral tree; every tree 

is either central or bicentral, but not both. 


Section 3 Greedy algorithms 


1 Spanning trees 

Let G be a connected graph. Ther. a spanning tree in G is 
a subgraph ot G that includes every vertex and is also a 
tree. 

Given a connected graph, a spanning tree mav be 
constructed by using either of the following methods. 
Building-up method 

Select edges of the graph one at a time, in such a way 


that no cycles are created; repeat this procedure until all 
vertices are included. 


Cutting-down method 
Choose any cycle and remove any one of its edges; 
repeat this procedure until no cycles are left. 


2 Minimum connector problem 

Let T be a spanning tree of minimum total weight in a 
connected weighted graph G. Then T is a minimum 
spanning tree or a minimum connector in G. 


Kruskal’s greedy algorithm for a minimum connector 
To construct a minimum spanning tree in a connected 
weighted graph G, successively choose edges of G of 
minimum weight in such a way that no cycles are 
created, until a spanning tree is obtained. 


Prim’s greedy algorithm for a minimum connecter 

To construct a minimum spanning tree T in a connected 

weighted graph G, build up T step by step as follows: 

¢ put an arbitrary vertex from the graph G into I; 

* successively add edges of minimum weight joining a 
vertex already in T to a vertex notin J, until a 
spanning tree is obtained. 

(Note that with Prim’s algorithm a connected graph is 

obtained at each stage.) 


Tabular (matrix) form of Prim's greedy algorithm for a 

minimum connector 

START with a table of weights for a connected weighted 
graph, and with no circled entries in the table 


SrePp 1 Delete all entries in column i, and mark row 1 
with 2 Star.” 


Step 2 Select a smallest entry from the uncircled entries 
in the row(s) marked with a star. 
If no such entry exists, STOP: the edges 
corresponding to the circled weights form a 
minimum connector, and its total weight is the 
sum of the circled weights. 
Otherwise, go to Step 3. 


STEP 3 (a) circle the weight w;j identified in Step 2; 
(b) mark row] witha se 
(c) delete the remaining entries in column j. 
Return to Step 2. 


| Theorem 3.1 


Prim's and Kruskal’s algorithms always produce a 
| spanning tree of minimum weight. 


3 Maximum connector problem 

Let T be a spanning tree of maximum total weight in a 
connected weighted graph G. Then T is a maximum 
spanning tree or a maximum connector in G. 


Prim’s greedy algorithm for a maximum connector 

To construct a maximum spanning tree T in a connected 

a graph G, build up T step by step as follows: 
¢ put an arbitrary vertex from the graph G into T; 

° successively add edges of maximum weight joining a 
vertex already in T to a vertex not in T, until a 
spanning tree is obtained. 


Kruskal’s greedy algorithm for a maximum connector 


Kruskal’s algorithm for a minimum connector can be 
adapted to give a maximum connector. 


4 Travelling salesman problem 
Given a weighted complete graph, find a minimum- 
weight Hamiltonian cycle. 


There is no known efficient algorithm for the travelling 
salesman problem — it is an NP-complete problem. 


One method for finding an approximate solution to the 
travelling salesman problem is to find a lower bound tor 
the total weight of a minimum-weight Hamiltonian cycle 
by solving a related minimum connector problem 
instead. 

Method for finding a lower bound for the solution to 
the travelling salesman problem 

STEP 1 Choose a vertex V and remove it from the graph. 
STEP 2 Find a minimum spanning tree connecting the 
remaining vertices, and calculate its total 

weight W. 

Find the two smallest weights, w, and wz, of the 
edges incident to V. 

STEP 4 Calculate the lower bound W + w, + w». 

Different choices of V give different lower bounds. The best 


lower bound is the largest, because it gives the most 
information about the actual solution. 


STEP 3 


Section 4 Multi-terminal flows 


1 The unit considers the following type of multi- 
terminal flow problem: 

find a maximum flow from any location in a given 
network to any other in the network, assuming that only a 
single pair of locations can communicate at any one 
time. 


2 Let N be a given undirected network. The complete 
network C(N) is the undirected network constructed by 
taking the complete graph with the same vertices as N, 
and assigning to each edge VW a capacity equal to f(VW), 
the value of a maximum flow between V and W (or W 
and V). 


; Theorem 4.1 | 
Let N be an undirected network, and let C(N) be the 

| corresponding complete network. Then, in any cycle 

| of C(N), the smallest capacity occurs more than once. 


3 A maximum spanning tree in C(N) can be used to 
reconstruct the capacities in C(N), using the capacity 
rule. 


Capacity rule 

To find f(VW) for two vertices V and W in an undirected 
network N, look at the path joining the vertices V and W 
in a maximum spanning tree of the complete network 
C(N); then f(VW) is the smallest capacity along this path. 


Theorem 4.2 | 


Let N be an undirected network with n vertices, and 
let C(N) be the corresponding complete network. 
Then C(N N) has at most n — 1 different capacities. 


single vertex. 


4 Two undirected networks are flow-equivalent if 

(a) they have the same number of vertices; 

(b) the maximum flows between pairs of vertices in the 
two networks are the same. 

A flow-equivalent tree is a cut tree if each branch of the 

tree corresponds to a minimum cut of the original 

network; that is, each branch, when removed, separates 

the vertices into two sets X and Y corresponding to a 

minimum cut in the original network. 


5 Algorithm of Gomory and Hu 

The algorithm is used for finding the maximum flows 
between all pairs of vertices in a given undirected 
network. 

The algorithm constructs a cut tree, branch by branch, by 
calculating a maximum flow for each branch. To attach 
each new branch to the existing tree, the cut tree 
property is used: the minimum cuts in the cut tree and the 
corresponding minimum cuts in the original network separate 
the same vertices. If the original network has n vertices, 
then a cut tree has n —1 edges, so only n-1 maximum 
flow calculations are required. 

Some of these calculations may be simplified by 
condensing certain vertices of the given network into a 


—— - iene, 


Condensing theorem 


Consider a minimum cut in an undirected network N 
separating the vertices into two disjoint sets X and Y. 


When finding a maximum flow between any two ' 
vertices in X, we may condense all the vertices in Y to 
a single vertex (and vice versa). ! 

' 


Edges connecting a vertex x in X to any vertex in Y are | 


_ replaced by a single edge with capacity equal to the ' 


sum of the capacities of the connecting edges. 


| YiyoUays | 


condensed 
vertex 


Once a cut tree has been constructed, the remaining 
maximum flows may be found by using the capacity rule, 
since a cut tree is a maximum spanning tree for the 
corresponding complete network. 
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Networks 2 Optimal paths 


Section 1 Algorithms using 
adjacency matrices 


1 


| Theorem 1.1 

Let Dbea digraph with n vertices labelled 1, 2, ..., 1; let 
A be its adjacency matrix with respect to this listing of 
_ the vertices, and let k be any positive integer. Then the 
number of walks of length k from vertex 7 to vertex j is 
the entry in row i and column j of the matrix Ak, 


Theorem 1.2 
Let D be a digraph with n vertices labelled 1, 2, ..., n; let 
- A be its adjacency matrix with respect to this listing of 
the vertices, and let B be the matrix 


B=At+Att+e+A7), 
Then D is strongly connected if and only if each non- 
diagonal entry in B is positive; that is, bj > 0 whenever 
LF}. 


— 
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2 Eulerian trails in connected digraphs 
An Eulerian trail in a connected digraph is a closed trail 
that includes every arc just once. 


Theorem 1.3 
_ A connected digraph is Eulerian if and only if, for each 
| vertex, the out-degree equals the in-degree. 


r 

Theorem 1.4 
_ An Eulerian digraph can be split into cycles, no two of 
| which have an arc in common. 


Algorithm for finding an Eulerian trail in a digraph 
To find an Eulerian trail in a digraph with no loops, and 
adjacency matrix A, carry out the following steps. 
START Set cycle counter p =0. 
STEP 1 Increase p by 1. 
Find a row k with a non-zero entry. 
Set n = k and put vertex vx at the start of a cycle: 
Cy = Uk... 
STEP 2 Choose m such that a, >0. 
Include vertex v,, as the next vertex in cycle C,,. 
Reduce 4,,,, by 1. 
Set 1 = m. 
If n xk, return to beginning of Step 2. 


If n= k and not all elements of A have been 
reduced to zero, store cycle C, and return to 
Step 1; otherwise go to Step 3. 
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STEP 3. If there is only one stored cycle, this is the 
required Eulerian trail: STOP. 


Otherwise, in the set of stored cycles, find two 
cycles C; and C; with at least one vertex v in 
common. In C;, replace the vertex v by C;,, after 
writing C; in the form C; = v...v. 


Delete C; from the set of stored cycles. 


Repeat this process until a single closed trail 
remains — this is the required Eulerian trail: 
STOP. 


3 Hamiltonian cycles in connected digraphs 
A Hamiltonian cycle in a connected digraph is a cycle 
which passes through every vertex. 
The algorithm for finding Hamiltonian cycles generates 
matrices in which each entry is a path or a cycle. The 
entries are called strings and the matrices are combined 
by latin multiplication (denoted by #), which is similar to 
matrix multiplication with the strings combined as 
follows: 
¢ string multiplication concatenates the strings; for 
example, 
VV 94 X U3V5 =U {UW 43s 
e string additions are written one above the other; for 
example, 
V0204 
U305U 6 
¢ a string which includes a vertex more than once is 0 
unless such a vertex occurs just at the beginning and 
the end; for example, 
V1V 90103 =0 
V1V7030, 40 
e 0 x (anything) = 0 
¢ 0 + (something) = (something) 


V0 704 i 03050 6 = 


Algorithm for finding Hamiltonian cycles in a digraph 
To find all the Hamiltonian cycles in a digraph with 
n vertices and no loops, carry out the following steps. 
STEP 1 Define an n Xn matrix C as follows. 
Put 
cj = vjv;, if there is an arc from 4; to ¥;, 
cj = 0, otherwise. 
Define an n x n matrix D to be the matrix 


obtained from C by deleting the first vertex in 
each non-zero entry of C. 


Set k = 1, where k is the power to which the matrix 


C is raised. 
STEP 2 Form C*t!= C*# D, where # denotes latin 
multiplication. 


STEP3 If k+1#n, increase k by 1 and return to Step 2. 
i k+le=7 SI1CP. 


The entries in C“ give the paths of length k, fork=1,..., 1. 
The entries in C’ give the Hamiltonian cycles. 


Section 2 Optimal path algorithms 


1 Shortest paths in weighted digraphs 

A shortest path from a vertex S to a vertex T ina 
weighted digraph is a path from S to T of smallest 
possible length. 


Shortest path algorithm 
This algorithm looks out from vertices of known potential. 
START Assign potential 0 to S. 
GENERAL STEP Consider the vertex (or vertices) just 
assigned a potential. 
For each such vertex V, consider each vertex W which 
can be reached from V along an arc VW, and assign W 
the label | 

(potential of V) + (distance VW) 
unless W already has a smaller label assigned from an 
earlier iteration. 
When all such vertices W have been labelled, choose the 
smallest vertex label in the network which is not already 
a potential, and make it a potential at each vertex where 
it occurs. 
Repeat the general step with the new potential(s). 
_ STOP when T has been assigned a potential. 
The shortest distance from S to T is the potential of T. 
To find a shortest path, trace backwards from T and 
include an arc VW wherever 

(potential of W) - (potential of V) = distance VW 
until S is reached. 


2 Longest paths in weighted digraphs 
A longest path from a vertex S to a vertex T in a weighted 
digraph is a path from S to T of greatest possible length. 


Longest path algorithm 
This algorithm looks back to vertices of known potential. 
START Assign potential 0 to S. 
GENERAL STEP Consider all vertices which can be 
reached directly only from vertices with known 
potentials. 
For each such vertex W, consider each arc VW into W, 
and assign W the label 

(potential of V) + (distance VW) 
unless W already has a larger label. 
When all such arcs VW have been considered, make the 
vertex label at W a potential. 
Repeat the general step with the new potential(s). 
. STOP when T has been assigned a potential. 


The longest distance from S to T is the potential of T. To | 


find a longest path, work backwards from T and include 
an arc VW wherever 

(potential of W) — (potential of V) = distance VW. 
If the process does not terminate — for example, if the 
network has.a cycle — there is no longest path 


Section 3 Critical path analysis 


This section is concerned with scheduling a project 
involving several activities when there 1s no restriction on the 
number of workers available. 

The precedence relations show which activities must 
precede which. 


1 Activity networks 
An activity network is a network which can be used to — 
schedule a project. It represents the activities from START 
to FINISH in accordance with the precedence relations. 

In an activity network where vertices represent activities, . 
each activity is represented by a numbered vertex. An 
arc joining a vertex X to a vertex Y indicates that the 
activity represented by X must be completed before the 
activity represented by Y can be started. 


Activity network construction algorithm 
Part A: procedure for numbering the vertices 


START Represent each activity by a vertex. 
For each vertex, create a shadow vertex, so that 
for each activity there are two corresponding 
vertices — the original vertex and the shadow 
vertex. 
Construct a bipartite graph in which one set of 
vertices consists of the original vertices, and the 
other set consists of the shadow vertices. 
If an activity Y must follow an activity X, draw an 
edge joining the original vertex representing Y to 
the shadow vertex representing X. 


al . 

Y OY 

Step 1 Number consecutively all the original vertices 
(chosen in any order) which have no edges 
incident with them. Record the numbering, 
together with the iteration number. 

STEP 2 Delete all numbered vertices, their 
corresponding shadow vertices, and all edges 
incident with these vertices. 

If not all the vertices have been numbered, return 
to Step 1. 


If all the vertices have been numbered, go to 
Part B, Step 3. | 
Part B: procedure for drawing the activity network 
STEP 3 Draw a START vertex, and the vertices numbered 
in the first iteration. 
Draw an arc from the START vertex to each vertex 
which was numbered in the first iteration. 
Assign a weight of zero to each arc. 
STEP 4 Draw the vertices which were numbered in the 
next iteration. 
To each such vertex Y, draw an arc from each 
previously numbered vertex X if there is an edge 
joining the original vertex Y to the shadow 
vertex X in the original bipartite graph 
constructed in Part A. 
Assign a weight to each arc XY equal to the 
duration of the activity X. 
Repeat until all vertices have been included in 
the activity network. 


STEP 5 Draw a FINISH vertex. 


From each terminal vertex Z (that is, each vertex 
whose out-degree is zero), draw an arc to the 
FINISH vertex. 
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Assign a weight to each such arc equal to the 
duration of the corresponding activity Z. STOP. 
The activity network has now been constructed. 


An alternative type of activity network uses arcs to 
represent activities and vertices to represent events. An 
event is a stage in the process corresponding to the start 
or finish of one or more activities. In this type of activity 
network, it may be necessary to introduce a dummy 
activity to show that an activity cannot start until another 
activity has been completed. Dummy activities have no 
time duration. They are also used when two activities 
have the same start and end events (each activity can 
then be referred to by a different number). 


2 Critical paths 

The minimum completion time for a project can be 
calculated from the activity network for the project. To do 
this, we find a longest path, called a critical path, from 
the START vertex to the FINISH vertex. By a longest path we 
mean a path for which the sum of the times associated 
with the arcs of the path (that is, the length of the path) is 
the greatest possible. The minimum completion time is 
equal to the length of a critical path. 

The maximum time by which an activity may be delayed 
without delaying the project is called the float of that 
activity. 


| Conventions 

_ For anetwork involving n activities, the START vertex is 

regarded as the Oth vertex and the FINISH vertex as the 
(n + 1)th vertex. 


The duration of activity i represented by the arc jj is 
denoted by ¢; ;. 
Ci 

—"~e 

: ] 
The algorithm assigns labels p; and e; to each vertex j 
for j = 0, 1,.., 1 +1. When the algorithm has been 
_ completed: 


eis the length of the longest path from the START 
' vertex to vertex 7; 


p.is the number of the preceding vertex on this 
longest path. 


——oe aig 


Critical path construction algorithm 

Part A: labelling procedure 

STEP 1 Assign to the START vertex (vertex 0) the labels 
p, =O ande,=0. 

STEP 2 Carry out the following procedure for each 
vertex j, starting with j = 1 and continuing with 
j =2, 7 =3, and so on, until all the vertices 
(including the FINISH vertex, corresponding to 
j =n +1) have been labelled. 
For the current vertex j, calculate, for each arc ij 
incident to vertex j, the sum 

ej; + Chip. 

Choose the maximum value of these sums for all 
such arcs ij, and set ¢ equal to this value. 
Set p. equal to the value of i for which this sum is 
largest; in the case of a tie, choose any of the 
appropriate values. 

Part B: tracing back procedure 

“Step 3 Start with the FINISH vertex n + 1, and mark the arc 

joining this vertex to the preceding vertex p,,+1. 


Srep 4 Consider the vertex j from which the last marked 
arc is incident. Mark the arc joining this vertex to_ 
the preceding vertex pj. 


18 


Repeat until the START vertex is reached. STOP. 
The marked arcs form a critical path. 


The sum of the weights on the arcs of the critical path is 
the minimum completion time, and is given by the value 
of Cy +1: 
3. Earliest and latest starting times 
The critical path construction algorithm assigns a label e; 
to each vertex i. The value of e; is the length of the longest 
path from the START vertex to the vertex i. For each 
activity i, e; is the earliest starting time of activity 7, since 
activity 7 cannot be started until all preceding activities 
have been completed. For each activity i, the latest 
starting time /; is the latest time at which activity i can be 
started without delaying the whole project. 
The float of activity i (the maximum amount of time by 
which the activity can be delayed without delaying the 
project) is the difference between these two times: 

float of activity 1 =1;- @;. 
For an activity i on the critical path, e; =1;. The float of an 
activity on the critical path is therefore zero. 


Algorithm for calculating latest starting times 
This algorithm is applied to an activity network (with 
n activities) for which the minimum completion time has 
been calculated. It is used to calculate the value of the 
latest starting time |; for each vertex i, for 1=n+1,N, ..., 
1 
STEP 1 For vertex n + 1 (the FINISH vertex), set 1,,.; equal 
to the minimum completion time. 
STEP 2 Carry out the following procedure for each 
vertex i, starting with i = n and continuing with 
i=n—1,i=n-—2, and so on, until all the vertices 
including vertex 0 (the START vertex) have been 
considered. 
For the current vertex i, calculate, for each arc ij 
incident from vertex i, the difference . 
Lj — Cj jr 
where c; ; is the duration of activity 7. 
Choose the minimum value of these differences 
for all such arcs ij, and set 1; equal to this value. 
When the START vertex has been considered, 
SlOr. 
(Note that 1) must be zero — if this value is not 
obtained, a mistake has been made somewhere.) 


Section 4 Scheduling 


1 An ideal schedule satisfies the following conditions, 

called the factory rules. 

1 No worker may be idle if there is some activity which, 
can be done. 

2 Oncea worker starts an activity, that activity must be 
continued by the same worker until it is completed. 

3. The project must be completed as soon as possible 
with the manpower available. 


2 Algorithm for scheduling activities 

The algorithm assigns activities to a number of 
processors, which may be people or machines. It is 
assumed that the activity network has been constructed 
and the latest starting times calculated: 

The algorithm makes use of a hypothetical clock — 
called the project clock — to keep track of time. At the 
end of each iteration, the project clock is advanced so 
that it records the time for which the project has been 
running. 


Critical path scheduling algorithm 
START Set the project clock to 0. 


STEP 1 If at least one processor is free, assign to any free 
processor the most critical unassigned activity 
which can be started. (The most critical activity is 
one with the smallest latest starting time.) 


Repeat until no processor is free or until no 
unassigned activity can be started. 


Step 2 Advance the project clock until a time is reached 
when at least one activity has been completed, so 
that at least one processor is free. 


If not all the activities have been assigned, return 
to Step 1, 

If all the activities have been assigned, advance 
the project clock until all the activities have been 
completed. STOP: the project clock gives the 
minimum completion time. 


3 Outline of algorithm incorporating a protection 

scheme for type of example in Section 4.2 

Step 1 Compute for each activity the sum of the earliest 
starting time and the latest starting time, and 
rank the activities in ascending order according 
to the values of these sums. 


STEP 2 Assign activities successively to free workers 
according to this ranking, taking into account the 
precedence relations. If a preceding activity has 
not already been assigned, break the ranking 
order when assigning the activities. 


Section 5 Bin packing 


A problem to determine the minimum number of 
workers required to complete a project within a given 
period is an example of a bin-packing problem. 

A bin-packing problem involves trying to pack items 
(activities) into the minimum number of ‘bins’, each of 
the apeced jul 


I eas fit eaciang digonihii 
Always place the next item to be packed in the current | 
bin if possible; otherwise, place it in the next bin. : 


Hs aes a 


First-fit packing algorithm | 
Always place the next item to be packed in the lowest 
numbered bin which can accommodate that item. 


Lap eet rr a a an te 


First-fit decreasing packing algorithm 
Order the items in decreasing order of size, and apply 
the first-fit procedure to this reordered list. 


rn a te 


If the items to be packed are available only one at a time, 
so no reordering of items is possible, an algorithm such 
as the next-fit algorithm or the first-fit algorithm must be 
used. Such algorithms are called on-line algorithms. 

If it is possible to reorder the items to be packed, the 
first-fit decreasing algorithm can be used. Such an 
algorithm involving reordering is called an off-line 
algorithm. 
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Design 2 Kinematic design 


Section 1 Kinematic structure 


1 Akinematic system is a mechanical system designed 
to produce, transmit, control, constrain, or resist 
movement. 

The design of kinematic systems tries to achieve various 
degrees of movement. In particular, it concerns the 
arrangment and interconnection of the component parts of 
a system so that their motions are suitably controlled and 
constrained. This is what is meant by kinematic 
structure. | 

The primitive rigid components of kinematic systems are 
called links. The interconnections between links are 
called joints. 

. To emphasize the fact that binary joints involve just two 
links, they are referred to as kinematic pairs. 

The six Reuleaux pairs are illustrated below. 


ee es Ea St te Er AE 


pair sketch contact surface 


revolute pair 


prismatic pair 


elliptical cylinder 


screw pair 


helicoid 


circular cylinder 


planar pair 


plane 


spherical pair 


sphere 
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2 The direct graph of a kinematic system is a graph in 
which the vertices represent the joints of the system and 
the edges represent the links of the system. In a direct 
graph, two vertices are joined by an edge if the 
corresponding joints belong to the same link. 

The ‘interchange graph of a kinematic system is a graph 
in which the vertices represent the links of the system 
and the edges represent the joints of the system. In an 
interchange graph, two vertices are joined by an edge if 
the corresponding links are connected at a common 
joint. 

3 The combined position and orientation of a link in 
space, relative to a reference coordinate system, is called 
its pose. If the pose changes with time, then the link is in 
motion. 


4 Types of kinematic system 


When all the axes of.a system are parallel, we refer to the 
system as planar because the links all move in parallel 
planes. When all the axes intersect in a single point, the 
system is called a spherical system, because the links all 
move on the surfaces of concentric spheres. If the axes 
are all mutually skew, the system is called a 
spatial system, because all the links have eeneral 


motions in space. 


5 The mobility of a kinematic system is the number of 
independent quantities or coordinates that must be 
specified in order to describe the pose of every link of the 
system with respect to a coordinate frame fixed on a 
chosen reference link, called the fixed link. 


Section 2 Braced rectangular 
frameworks 


1 A rectangular framework consists of one or more 
rectangles, or bays. 

A brace is a mechanical restriction on the motion of a 
bay in a rectangular framework which reduces the bay’s 
mobility by 1, so that one fewer quantity is required to 
specify the pose of every link of the framework. 


2 Ina standard rectangular framework, a row is the set. 
of all links forming the vertical sides of a horizontal string 
of bays, and a column is the set of all links forming the 
horizontal sides of a vertical string of bays. 
If the rows are numbered 1, fr, ..., %, ..., % sequentially 
from top to bottom and the columns are numbered c,, ¢:, 
» Gy ++ Cn Sequentially from left to right, then bay (i, j) is 
the bay whose vertical links belong to row r; and whose 


horizontal links belong to column c;. 


Theorem 2.1 | 
A braced rectangular framework, with rows n, 1, ..., f, | 
ft, and columns €,, Cy... Cy -++ Cyn 18 rigid if and onlv if 
its braces are located such that, for i= 1,2,..., 1 and 

j=1,2,...,m 


(a) r; must remain parallel to 1; for all r, and r;; 


(b) c; must remain parallel to c;, for all c; and c; 
(c) 4+; must remain perpendicular to c;, for all r; and c;; 
under any attempted deformation of the framework. 


L i ps i i 
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2 Theorem 2. 2: | 


| 
| A braced rectangular framework is rigid if and only if 
| its associated bipartite graph is connected. | 


ee 


A bracing of a framework is a particular Niseation of 
braces to bays of the framework. 


4 Arigid braced rectangular framework is minimally 
braced if no brace can be removed without destroying 
the rigidity; the corresponding bracing is a minimum 
bracing. 


Theorem 2.3 | 


A rigid braced rectangular framework is minimally — 
braced if and only if its associated bipartite graph is a | 
spanning tree. | 


Corollary 

| If the bipartite graph associated with a rigid braced 
rectangular framework has either of the following 
properties, then the corresponding bracing is not a | 
minimum bracing: | 


e the graph has n vertices and more than n - 1 
edges; 
e the graph contains a cycle. 


ee eee 
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Section 3 Freedom and constraint 


1 The freedoms of a link are the independent 


quantities or coordinates that must be specified in order 
to describe the pose of the link with respect to some 
reference coordinate frame. 
The freedoms of a two-link kinematic system are the 
independent quantities or coordinates that must be 
specified in order to describe the pose of one link of the 
system with respect to some reference coordinate. frame 
fixed in the other link of the system. 
Then 3.1 
At least three points must be considered when | 
describing the freedoms of a link in space. | 
Furthermore, once the positions of any three distinct — 
_ non-collinear points of a link have been specified, — 
then the pose of the whole link has been ppeaies | 


2 Aconstraint on a link is a geometric restriction that 
removes one of its freedoms. The pattern is summarized 
by 

M =F-C, 
where M is the mobility, F is the number of freedoms 
and C is.the number of constraints. 


3 The connectivity of a kinematic pair considered as 
the only joint in a two-link kinematic system, is the 
number of freedoms of one of the links with respect to 
the other. 


Se pr anata a. sees 


Section 4 Planar kinematic systems 


In this section, attention is restricted to those planar 
kinematic systems ee ay coe pairs. 


— Ne Su. aera resase 


— 
| Theorem 4.1: first Hebe criterion 

_ Ina planar kinematic system containing only revolute 
pairs, the mobility M of the system relative to one of , 
the links (considered as the fixed link) is given by the | 
mobility criterion ; 


M =3(n-1)-3j 


where n is the total number of links (regardless of 
_ their multiplicity) and 7 is the number of revolute |, 
_ pairs. 3 
Pe 
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The value obtained for M from the mobility criterion in 
Theorem 4.1 is always an integer. There are the following 
three possible cases. 

If M is positive and if there are no geometrical 
complications, then the links of the planar system can 
have relative motion and the mobility of the joints must 
be controlled in order to control the motion of the 
system. The system is said to be mobile. 

If M is zero and if there are no geometrical 
complications, then the links of the system have no 
relative motion and the system is said to be immobile. 

If M is negative and if there are no geometrical 
complications, then the links of the system again have no 
relative motion. But, in addition, some of the links and 
joints can be removed without allowing the remaining 
links to have any relative motion. The system is said to be 
overconstrained. 


Theorem 4.2: second mobility criterion . 
In a planar kinematic system containing links and _ 
_ joints of various multiplicities, where each r-ary joint is : 
_ equivalent to r—1 revolute pairs, the mobility M of the 
system relative to one of the links (the fixed link) is 
given by the mobility criterion 


((2r —3)n, +++» +(2s—3)n,) -3, 


where g is the total number of joints (regardless of 
their multiplicity), no link has multiplicity greater than 
_ s,and n, is the number of r-ary links (2 <r<s)._ 


= eye ee eee ek are reer 


a ee eee ee aS 


| 
| 
M =2g+3n, +1, - 
| 


| 
| 


_ Theorem 4.3: mobility criterion for graphs 
In a simple connected graph whose vertices all have | 
degree greater than or equal to 2, and which therefore - 
can be considered to be the direct graph of a planar | 

kinematic system containing multiple joints (each of — 

- which can be regarded as a combination of revolute 
pairs) and only binary links, the mobility M of the 
graph relative to one of the edges (considered to be 
the fixed edge) is given by the mobility criterion 


M =2¢- 


_ where g is the total number of vertices of the graph | 
| (regardless of their degree) and n, is the number of | 
| edges of the graph. | 
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As with planar kinematic systems, the value obtained for 


M from the mobility criterion in Theorem 4.3 is always an 


integer. Again, M can be positive, zero or negative, and, 
by analogy with the planar kinematic systems, the graph 
is then said to be mobile, immobile or overconstrained, 
respectively. 


Graphs 3 Planarity and colouring 
Section 1 Planarity 


1 Planar graphs 

A graph G is planar if it can be drawn in the plane in 
such a way that no two edges meet except at a vertex with 
which they are both incident. Any such drawing is a 
plane drawing of G. 

If no plane drawing of G exists, then G is non-planar. 


K33;and K; are non-planar. 


2 Euler’s formula 
Let G be a planar graph. Then any plane drawing of G 
divides the set of points of the plane not lying on G into 
regions, called faces; one of these faces is of infinite 
extent and is called the infinite face. 

Let G be a planar graph, and let f be any face of a plane 
drawing of G. Then the degree of f, denoted by deg f, is 
the number of edges encountered in a walk around the 
boundary of the face f. 

If all faces have the same degree ¢g, then G is face- 
regular of degree g. 
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Theorem 1.1: handshaking lemma for planar graphs 
| In any plane drawing of a planar graph, the sum of all 
_ the face degrees is equal to twice the number of 
, edges. 
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' Theorem 1.2: Euler’s formula for planar graphs 
_ Let G be a connected planar graph, and let n, m and f 
_ denote, respectively, the numbers of vertices, edges 
and faces in a plane drawing of G. Then 


n—-mt+f=2, 
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_ Corollary 1.1 
_ Let G be a simple connected planar graph with n (2 3) 
| vertices and m edges. Then 


m <3n—-6. 
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Corollary 1.2 | 
' Let G be a simple connected planar graph with n (2 3) 
vertices, m edges and no triangles. Then | 


, m <2n-4. 
Corollary 1.3 | 
Let G be a simple connected planar graph. Then G 


contains a vertex of degree 5 or less. 


el 


Corollaries 1.1 and 1.2 can be used to prove (by 
contradiction) that certain graphs are non-planar. 


3 Subdivisions and contractions 
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insert vertices 
of degree 2 bg 


G subdivision of G 


Any graph formed from a graph G by inserting vertices 

of degree 2 into the edges of G is called a subdivision 

of G. 

¢ If Gisa planar graph, then every subdivision of G is 
planar. 

e If Gis a subdivision of a non-planar graph, then G is 
non-planar. 


22 


e If Gis a graph that contains a subdivision of K; or K33, 
then G is non-planar. 

e If Gis anon-planar graph, then G contains a 
subdivision of K; or K;3. 


nee: i 


em 


Theorem 1.3: Kuratowski’s theorem 


A graph is planar if and only if it does not contain a 
subdivision of K; or K33 ‘ 


Another characterization of planar graphs involves the 
notion of ‘contracting’ an edge. This is done by bringing 
one vertex closer and closer to the other vertex until they 
coincide, and then coalescing any resulting multiple 
edges into a single edge. 

A contraction of a graph is the result of a sequence of 
edge contractions. For example, K; is a contraction of the 
Petersen graph, since it is the result of contracting each 
of the five ‘spokes’. 


contract 
spokes 


G contraction of G 


Theorem 1.4 


A graph is planar if and only if it does not contain a 
subgraph that has K, or K33 as a contraction. 


Theorems 1.3 and 1.4 do not provide an easy way of 
showing that a given graph is planar, since this would 
involve looking at a large number of subgraphs and 
verifying that none of them is a subdivision of K; or K;;, 
or contains K, or K,; as a contraction. For this reason, no 
currently used algorithm for testing the planarity of a 
graph is based on either of these theorems. 


4 Duality 

(The following definition assumes that a plane drawing of 

G is available.) 

Let G be a connected planar graph. Then a dual graph 

G* is constructed from a given plane drawing of G, as 

follows. 

(a) Draw one new vertex 1n each face of the plane 
drawing — these are the vertices of G*. 

(b) For each edge e of the plane drawing, draw a line 
joining the vertices of G* in the faces on either side 


G G* 


Different plane drawings of a planar graph G may give 
rise to non-isomorphic dual graphs G*. 

If G is a plane drawing of a connected planar graph, then 
so is its dual G*, so (G*)*, the dual of G* can be 
constructed. 

The graph (G*)* is isomorphic to G. 
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Theorem 1.5 

Let G be a plane drawing of a connected planar graph 
with n vertices, m edges and f faces. Then G* has 
f vertices, m edges and n faces. 


| 


| 
| 
| 
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Correspondences are = given | in the following table. 
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- dual graph G* 


plane drawing Ge 


an edge of G 


an 1 edge c of G* 
a vertex of degree k in G a face of degree k in G* 
a face of degree k in G a vertex of degree k in G* 


a cutset of G* with k edges 
a cycle of length k in G* 


a cycle of length kin G 
a cutset of G with k edges 


These correspondences can be used to obtain new 
results from old ones by dualizing. 

Since loops (cycles of length 1) and pairs of multiple 
edges (cycles of length 2) in G correspond to cutsets with 
1 and 2 edges in G*, Corollary 1.1 can be dualized to give 
the following theorem. 


- Theorem 1.6 


| Let G* be a connected planar graph with f faces and 
_m edges, and with no cutsets with 1 or 2 edges. Then 


m S3f- 6. 
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Similarly, eoeien- 41 1.3 can be dualized to give the 


following theorem. 


| 
| 
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| Theorem 1.7 
Let G* be a connected planar graph with no cutsets _ 
| with 1 or 2 edges. Then G* has a face of degree 5 or | 
_ less. | 


5 Testing for planarity 

e If the graph is disconnected, then it is planar if and 
only if each of its components is planar. 

e If the graph has a cut vertex (a vertex whose removal 
disconnects the graph), then it is planar if and only if 
each of the subgraphs obtained when the graph is 
‘broken apart’ at the cut vertex is planar. 

e If the graph has loops or multiple edges, then it is 
planar if and only if the graph obtained by removing 
the loops and coalescing the multiple edges is planar. 


Cycle method for planarity testing 
Given a graph G to test for planarity, the idea is to find a 
Hamiltonian cycle, to draw this cycle as a _ regular 
polygon, and then to try to draw the remaining edges so 
that no edges cross. 
Choose a Hamiltonian cycle C. Then list the remaining 
edges of Gand try to divide them into two disjoint sets A 
and B, as tollows: 
A is a set of edges which can be drawn inside C without 
crossing; 
B is a set of edges which can be drawn outside C 
without crossing. 
If this is possible, the graph G is planar, and the sets A 
and B can be used to obtain a plane drawing of G. If this 
is not possible, the graph G is non-planar. 
Two edges are incompatible if they cannot both be 
drawn inside C, or both be drawn outside C, without 
crossing. Otherwise they are compatible. 


Section 2 Colouring maps 
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| Theorem 2.1: four colour theorem for maps 


_ The countries (faces) of any map can be coloured with | 
four (or fewer) colours in such a_ way that 
| _ neighbouring: countries are coloured differently. | 
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A map is a plane ene of a connected planar graph 
containing no cutsets with 1 or 2 edges (in particular, a 
map contains no vertices of degree 1 or 2). A country is a 
face of a map. 

Every map must contain at least one of the following: 

* acountry with just two boundaries — a digon; 

¢ acountry with just three boundaries — a triangle; 


¢ acountry with just four boundaries — a square; 

¢ acountry with five boundaries — a pentagon; 

¢ two adjacent pentagons; 

¢ a pentagon adjacent to a hexagon. 

Such a set of configurations is an unavoidable set. 

A configuration is reducible if, when it is removed, the 
remaining map four-coloured, and then the 
configuration reinstated, the colouring can always be 
siacticonas to that configuration. 


| _ Theorem 2.2: six colour ihewein for maps 


_ The countries of any map can be coloured with six (or 


| fewer) colours in such a way that neighbouring 


~ countries are re coloured differently. 
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Theorem 2.3: five colour theorem for maps 


_ The countries of any map can be coloured with five (or 
fewer) colours in such a way that neighbouring 
_ countries are coloured differently. 


Section 3 Vertex colourings and 
decompositions 


Attention is restricted to simple graphs. 


1 Vertex colourings 

Let G be a simple graph. A k-colouring of G is an 
assignment of k colours to the vertices of Gin such a way 
that adjacent vertices are assigned different colours. If G 
has a k-colouring, then G is k-colourable. 

The chromatic number of G, denoted by x(G) is the 
mae number k Or which Gi is 3 k- colourable. 
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| phere Fe | 
Let G be a simple graph whose maximum _ vertex 
_ degree is d. Then 


¥(G) <d+1. 


Theorem 3.2: Brooks’ theorem : 
_ Let G be a connected simple graph whose maximum 
vertex degree is d. If Gis neither a cycle graph withan » 
odd number of vertices, nor a complete graph, then 
KG): = d. 


To find the chromatic number rxlG) 0 of a seanvle ah G 

Try to find an upper bound and a lower bound which are 

the same; then y(G) is equal to this common _ value. 

(If ¥(G) <k and y(G) 2k, then x(G) = k.) 

a upper bounds for x(G) 
the number of colours in an explicit vertex colouring 
of G 

e the number 7 of vertices in G; 

e d+1,where d is the maximum vertex degree in G: 

e d, where d is the maximum vertex degree in G, 
provided that Gis not C,, (for odd n) or K,,. 

possible lower bound for x(G) 

¢ the number of vertices in the largest complete 
subgraph in G. 

2 Colouring planar graphs 
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| | Theorem 3.3: six colour theorem for iewar araphs 


| 

The vertices of any simple connected planar graph 
can be coloured with six (or fewer) colours in such a | 
| way that adjacent vertices are coloured differently. | 


| | Theorem 3.4: dive volar ieoren for planar graphs 


The vertices of any simple connected planar graph 
can be coloured with five (or fewer) colours in such a 
| way that adjacent vertices are coloured differently. 
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| Theorem 3.5: four colour theorem for planar graphs 


The vertices of any simple connected planar graph 
can be coloured with four (or fewer) colours in such a 
way that adjacent vertices are coloured differently. 


3 Finding a vertex colouring 
Greedy algorithm for vertex colouring 
START witha graph G and list of colours 1, 2, 3, ... 


STEP 1 Label the vertices a, b,c, ...in any manner. 


Step 2 Colour the vertex labelled with the earliest letter 
in the alphabet not yet coloured with the first 
colour in the list not used for any adjacent vertex 
already coloured. 


Repeat Step 2 until all the vertices are coloured, then 
STOP. 

A vertex colouring of G has been obtained. The number 
of colours used depends on the labelling chosen for the 
vertices in Step 1. 


: iieoeent 3.6 

| For any graph G, there is a labelling of the vertices for 
| which the greedy algorithm yields a vertex colouring 
| with x(G) colours. | 


| 
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4 Vertex decompositions 

A vertex decomposition of G is obtained by splitting the 
set of vertices of G into disjoint subsets. 

A dominating set of vertices in a graph G is a set S of 
vertices with the property that each vertex of G is either 
in S or adjacent to a vertex of S. 

A minimum dominating set is a dominating set of 
smallest possible size. 

The dominating number of G, denoted by dom (G), is 
the number of vertices in a minimum dominating set. 


For colouring problems, in each subset, no two vertices are 
adjacent. For domination problems, each subset contains a 
vertex adjacent to all the other vertices in that subset. 

The independence problem for a graph G is that of finding 
the largest possible set of vertices of G, no two of which 
are adjacent. 


An independent set of vertices in a graph G is a set of 
vertices of G, no two of which are adjacent. 

A maximum independent set is an independent set of 
largest possible size. 

The independence number of G, denoted by ind (G), is 
the number of vertices in a maximum independent set. 
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| Theorem 3.7 


: For any graph G with n vertices: 
_ (a) dom (G) Sind (G); 
_ (b) x(G) x ind (G) 2 1. 


Section 4 Edge colourings and 
decompositions 


Attention is restricted to graphs without loops. 


1 Edge colourings 

Let G be a graph without loops. A k-edge colouring of G 
is an assignment of k colours to the edges of G in such a 
way that any two edges meeting at a vertex are assigned 
different colours. If G has a k-edge colouring, then G is k- 
edge colourable. 

The chromatic index of G, denoted by y’(G), is the 
smallest number k for which G is k-edge colourable. 


| Theorem 4.1: Vizing’s theorem 
Let G be a simple graph whose maximum vertex | 
_ degree is d. Then 

d<x(G)<d+l. 
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Theorem 4.2: Vizing’s theorem (extended version) 


Let G be a graph whose maximum vertex degree is d, 
and let h be the maximum number of edges joining a 
pair of vertices. Then 


d<¥/(G)<d+h. 
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Theorem 4.3: Shannon’s theorem 


| Let G be a graph whose maximum vertex degree is d. 
| Then 


dsx(G)s 
d=x(G)s 


‘ 


SaiZ, if d is even; 
(3d-1)/2, if dis odd. 
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To find the chromatic index y’(G) of a graph G without 

loops 

Try to find an upper bound and a lower bound which are 

the same; then y’(G) is equal to this common value. 

(If ¥’(G) <k and y’(G) 2k, then y’(G) = k.) 

aaa upper bounds for y’(G) 
the number of colours in an explicit edge colouring 
of G; 

e the number m of edges in G; 

e d+1, where d is the maximum vertex degree in G, 
provided that G has no multiple edges; 

e d+h, where d is the maximum vertex degree in G and 
h is the maximum number of edges joining a pair of 
vertices; 

e 3d/2, where d is the maximum vertex degree and d is 
even; 

e (3d-—1)/2, where d is the maximum vertex degree and 
dis odd. 

possible lower bound for y’(G) 

e d, the maximum vertex degree in G. 


2 Classifying graphs 


Theorem 4.4 
For the complete graph K,, | 


¢(Kn)= Qo 


n if nis odd. | 


if 11s even; 


Theorem 4.5: Kénig’s theorem 


Let G be a bipartite graph whose maximum vertex | 
degree is d. Then | 
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3 Finding an edge colouring 

Greedy algorithm for edge colouring 

START with a graph G and list of colours 1, 2, 3, .... 

STEP 1 Label the edges a, b, c, ... in any manner. 

STEP 2 Colour the edge labelled with the earliest letter in 
the alphabet not yet coloured with the first colour 


in the list not used for any previously coloured 
edge that meets it at a vertex. 


Repeat Step 2 until all the edges are coloured, then 
STOP. 

An edge colouring of G has been obtained. The number 
of colours used depends on the labelling chosen for the 
edges in Step 1. 


Theorem 4.6 
For any graph G, there is a labelling of the edges for 
_ which the greedy algorithm yields an edge colouring | 


| _ with x"(G ) colours. 


4 Edge decompositions 
An edge decomposition of G is obtained by splitting the 
set of edges into disjoint subsets. 

A matching in a graph G is a set of edges of G, no two of 
which have a vertex in common. 

The thickness of a graph G, denoted by f(G), is the 
minimum number of planar graphs that can be 
superimposed to form G. 


Notation 
Let a be any positive number. Then La is the integer 
obtained by ‘rounding a down’, and| a] is the integer 
obtained by ‘rounding 4 up’. 

The connection between these functions is given by 


a/b |=L(a/b) + (b-1)/ 0). 
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Theorem 4.7 


Let G be a simple connected graph with n (2 3) 
vertices and m edges. Then 


(a) t((G)=/m/(Gn-6)|; 
 (b) t(G)2 [m /(2n—4)]|, if G has no triangles. 


| Theorem 4.8 
(a) HK,)>L(n+7)/6]; 
(b) t(K,;) 2] rs/(2r+2s—4)]. 
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5 Decomposition into spanning subgraphs 


Some problems reduce to that of decomposing a given 
graph G into the maximum number of connected 
subgraphs, each of which includes every vertex of the 
graph. This number is denoted by s(G). 
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| Theorem 4.9 
| Let G be a connected graph with n vertices. Then s(G) 
is the largest integer for which the following statement 
is true: 
for each positive integer k <n, at least (k—1) x s(G) 
edges must be removed in order to disconnect G 
into k components. 
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Theorem 4.10 
Let G be a connected graph with n vertices and s(n — 1) 
edges. Then G can be decomposed into s spanning 
| trees if and only if 
for each positive integer k <n, at least (k-—1) xs 
edges must be removed in order to disconnect G 
into k components. 
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Networks 3 Assignment and 
transportation 


1 Matching problem 


1 Matchings 

The pairing of some or all of the elements of one set with 
elements of a second set is called a matching. In a 
matching in a graph, the elements are the vertices, and a 
matching is indicated by thick edges. 

A maximum matching is a matching with the largest 
possible number of edges. 


2 Marriage problem 

Given a set of men, each of whom knows some women 
from a given set of women, under what conditions is it 
possible for all the men to marry women they know? 


Theorem 1.1: marriage theorem | 
A necessary and sufficient condition for there to be a 
solution to the marriage problem (that is, for every 
man to be able to marry a woman he knows) is: 
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| for every subset of m men, the m men collectively _ 
know at least m women, for all values of m in the ! 
range 1 < m <n, where n is the total number of 
men. 


marriage condition. 


3 Modified marriage problem 
When it is not possible for all the men to marry, it is 
natural to ask the following question, called the modified 
marriage problem: 
what is the maximum number of men who can marry 
women they know? 


If the marriage condition is not satisfied, it can be written 

instead that n — d men can marry women they know if 
(number of women known by each subset of m men) 2 
m —d, for all values of m in the range 1 < m <n, and for 
some positive integer d. 

This is called the modified marriage condition. 


Theorem 1.2: modified marriage theorem 
If a group of nm men each know some of a group of | 
women, the maximum number of men who can marry | 
women they know is equal to the minimum value of 
the expression: 


(number of women known by a subset of m men) + _ 
(n— m), for any subset of m men, and for all values 
of m inthe range 1<m<n. 


2 Maximum matching problem 


1 Alternating paths 

Let G be a bipartite graph in which the set of vertices is 

divided into two disjoint subsets X and Y. An alternating 

path with respect to a matching M in G is a path which 

satisfies the conditions: 

(a) the path joins a vertex x in X to a vertex yin Y; 

(b) the initial and final vertices x and y are not incident 
with an edge in M; 

(c) alternate edges of the path are in M, and the other 
edges are notin M. 
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If an alternating path with respect to a matching M can 
be found, then a matching M’ which has one more edge 
than M can be constructed as follows. 


To find a new matching with one more edge 

Find an alternating path P with respect to an existing 
matching M. 

Form a new matching M’ with the following edges: 


e the edges of the alternating path P not in the 
matching M 


and 


e the edges of the matching M not in the alternating 
path P. 7 


2 Maximum matching problem 
How can a maximum matching in a bipartite graph be 
found? 


The maximum matching algorithm is given opposite. 


3 Assignment problem 


1 Assignment problem 

If a bipartite graph represents a practical problem, some 
of the matchings may represent more desirable solutions 
to the problem than others. This fact is expressed by 
assigning costs to the edges of the bipartite graph. These 
costs are usually displayed in the form of a cost matrix. 
How can a maximum matching with the lowest total cost 
be found? 

The Hungarian algorithm for the assignment problem 
involves modification of the cost matrix. For a given cost 
matrix, the graph containing the original vertices, 
together with the set of edges whose current cost is zero, 
is called the corresponding partial graph. 

The Hungarian algorithm for the assignment problem is 
given on page 28. 


Use of dummy vertices 

The Hungarian algorithm for the assignment problem 
applies to a bipartite graph with two sets of vertices X 
and Y containing the same number of elements. Given a 
problem in which this condition is not satisfied, it is 
possible to solve the problem by including dummy 
vertices and corresponding high entries in the cost 
matrix. The algorithm can then be applied in the usual 
way. 


2 Bottleneck assignment problem 

Suppose that a product is to be made on a serial 
production line which involves a number of activities, 
each of which must be completed once to produce one 
item of the product. A number of people are to be 
assigned to these activities. It is known how long it takes 
each person to complete each activity. How can the 
people be assigned to the activities so that the time 
taken to complete one item of the product is as short as 
possible? 

Algorithm for bottleneck assignment problem 

START with a given cost matrix. 


set p = (0. 


STEP 1 (a) Set d = smallest non-zero entry in the current 
cost matrix. 


(b) Reduce each non-zero entry in the current 
cost matrix by d. 


(c) Increase p by d. 
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Step 2 Use the matching algorithm to find a maximum 
matching in the bipartite graph. 


If a complete assignment is achieved, STOP: this 
is an optimum assignment and the value of p is 
the required answer. 


Otherwise, return to STEP 1. 


In simple examples a maximum matching is usually 
found by inspection. In more complicated examples, in 
Step 2 it is necessary to find an alternating path with 
respect to the existing matching and use it to obtain a 
maximum matching. 


4 Transportation problem 


1 Transportation problem 

A manufacturer has a number of factories, each of which 
can supply a particular product to a number of 
warehouses. Each factory has a fixed output, and there is 
a transportation cost involved in sending the product 
from each factory to each warehouse. To which 
warehouses should the products of each factory be sent 
so that the requirements of all the warehouses are 
satisfied at minimum total cost? 

In the bipartite graph for this type of problem, the 
factories are represented by black vertices, called supply 
vertices. The warehouses are represented by white 
vertices, called demand vertices. 

The Hungarian algorithm for the transportation problem 
is given on page 30. 

Use of a dummy vertex 

In a transportation problem in which the total supply is 
greater than the total demand, it is possible to use the 
Hungarian algorithm by introducing a dummy vertex. 
This represents an artificial warehouse (or other demand 
vertex) whose demand is chosen to be equal to the 
excess supply. 

In the cost matrix, the transportation costs from each 
supply vertex to this dummy vertex must be equal. After 
the algorithm has been used, a flow pattern which is a 
solution to the original problem can be found by 
removing all flows to the dummy vertex. 


2 Transhipment problem 

In the transportation problem it is assumed that the 
products from each factory are sent directly to one or 
more warehouses. However, in a practical situation, it 
may be cheaper to send the products from a factory first 
to an intermediate point (one of the other factories or 
warehouses) and then on to the destination warehouse. 
The sending of goods via intermediate points is called 
transhipment, and a problem in which transhipment is 
allowed is called a transportation problem with transhitpment, 
or, more simply, a transhipment problem. The graphical 
representation of a transhipment problem is not a 
bipartite graph, since goods can be sent from factory to 
factory, or from warehouse to warehouse. However, the 
graphical representation can be transformed into a 
bipartite graph by representing the factory or warehouse 
by two vertices, one representing the factory or 
warehouse as a supply point, and one representing it as a 
demand point. 

It is necessary to add to all the demands and supplies an 
amount of demand or supply at least as great as the total 
demand or supply. The algorithm for the transportation 
problem is then applied in the usual way. When the 
procedure of the algorithm has been completed, it is 
necessary to remove the artificial demands and supplies 
(and the corresponding artificial flows) to obtain a 
solution to the original problem. 


Algorithm for finding a maximum matching 


START with any matching M in a bipartite graph whose sets of vertices 
are X = {X1, Xo, ..., X,} and Y = {yj, Yo, .-., Ym}- 


Part A: labelling procedure 
STEP 1. Label with (*) each vertex in X which is not incident with any edge 
in the current matching M. 


If no vertex in X can be labelled, STOP: the current matching is a 
maximum matching. 


Otherwise, go to Step 2. 


STEP 2. Choose a newly labelled vertex in X, say x; and label with x; all 
the unlabelled vertices in Y joined to x; by an edge NOT IN M. 


Repeat this for all newly labelled vertices in X, and then go to 
Slep 3. 


STEP 3 Choose a newly labelled vertex in Y, say y; and label with (yj) all 
unlabelled vertices in X joined to y; by an edge IN M. 


Repeat this for all newly labelled vertices in Y. 
Repeat Steps 2 and 3 until 


EITHER a vertex in Y which is not incident with an edge in M is labelled 
(this is called breakthrough), in which case go to Part B, Step 4, 


OR no such vertex exists, in which case STOP: the current matching is 
a maximum matching. 


Part B: matching improvement procedure 


STEP 4 Find an alternating path as follows. 


Start at the breakthrough vertex, and go to the vertex indicated by 
its label. From this vertex, go to the vertex indicated by its label, 
and so on, until a vertex labelled (*) is reached. 


This path P is an alternating path. 
STEP 5 Form anew matching from: 


e the edges of the current matching M which are NOT IN the 
alternating path P; 


e the edges of the alternating path P which are NOT IN the 
current matching M. 


Remove all labels and return to Part A, Step 1 to see whether the 
new matching can be improved further. 


Summary of algorithm 
START with any matching. 


Part A: labelling procedure 


Label the vertices to identify an alternating path. 
If breakthrough is achieved, go to Part B. 


If breakthrough is not achieved, STOP: the current matching is a maximum 
matching. 


Part B: matching improvement procedure 
Find an alternating path by tracing back through the labels. 


Form a new matching from: 

e the edges in the current matching NOT IN the alternating path, 
e the edges in the alternating path NOT IN the current matching. 
Return to Part A. 


There is no difficulty in finding an 


initial matching M. If necessary, it can 


consist of just one edge. 
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Hungarian algorithm for the assignment problem 


START with no initial matching in a bipartite graph whose sets of vertices 
are X = {X1,X2,...,X,} and Y= {Y1,Y2, Sey Yn}. 


Construction of first partial graph 


STEP 0a Assign weights to the vertices as follows. 


To each vertex of X, assign a weight equal to the lowest cost on 
any edge incident to that vertex. 


Decrease the cost on each edge by the weight on the vertex to 
which it is incident. 


Repeat this procedure for each vertex in Y. 


STEP Ob Construct a partial graph consisting of the vertices of the original 
bipartite graph, together with only those edges whose current cost 
iS Zero. 


Part A: labelling procedure 
STEP 1 Label with (*) each vertex in X which is not incident with any edge 
in the current matching M. 


If no such vertex exists, STOP: the current assignment is an 
optimum assignment. 


Otherwise, go to Step 2. 


STEP 2. Choose a newly labelled vertex in X, say x; and label with (x,) all 
the unlabelled vertices in Y joined to x; by an edge NOT IN M. 


Repeat this for all newly labelled vertices in X. 


If no vertex in Y is labelled, STOP: the current assignment is an 
optimum assignment. 


Otherwise, go to Step 3. 


STEP 3. Choose a newly labelled vertex in Y, say y; and label with (y;) all 
unlabelled vertices in X joined to y; by an edge IN M. 


Repeat this for all newly labelled vertices in Y. 
Repeat Steps 2 and 3 until 
EITHER a vertex in Y which is not incident with an edge in M is labelled 
(this is called breakthrough), in which case go to Part B, Step 4, 
OR no more labelling is possible, in which case go to Part C, Step 6. 


Part B: matching improvement procedure 


STEP 4 Find an alternating path as follows. 


Start at the breakthrough vertex, and go to the vertex indicated by 
its label. From this vertex, go to the vertex indicated by its label, 
and so on, until a vertex labelled (*) is reached. 


This path P is an alternating path. 
STEP 5 Form anew matching from: 


e the edges of the current matching M which are NOT IN 
the alternating path P; 


e the edges of the alternating path P which are NOT IN the 
current matching M. 


Remove all labels and return to Part A, Step 1 to see whether the 
new matching can be improved further. 


Part C: modification of partial graph procedure 

STEP 6 Find the lowest current cost on any edge of the original bipartite 
graph which 
e starts ata LABELLED vertex in X; 
e ends at an UNLABELLED vertex in Y; 
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An original edge cost can be thought of 
as being the total cost of a trip from 
one vertex to the other. This cost is 
divided into three parts: the cost of 
leaving the vertex in X; the cost of the 
journey (the new edge cost); and the 
cost of arriving at the vertex in Y. 


When Step 0a has been carried out, for 
any edge, the sum of the weights on 
the end vertices plus the new cost is 
equal to the original cost. 


Note that the first time Step 1 is met, 
the current matching has no edges. 


e has anon-zero current cost. 
Call this lowest cost 6. 


STEP 7 (a) Increase the weight on each labelled vertex in X by 6. Note that in Step 7 no negative costs 
are produced in the cost matrix, since 6 
was chosen as the lowest non-zero 


(b) Decrease the weight on each labelled vertex in Y by 6. 


(c) For each edge of the original bipartite graph which joins a cost. Also, the total cost of a ‘trip’ from 
labelled vertex in X to an unlabelled vertex in Y, decrease the —_a vertex in X to a vertex in Y via any 
current cost by 6. edge not in the partial graph, that is, 


the sum of the two vertex weights and 
the edge cost, is kept constant 
throughout Step 7. 


(d) For each edge of the original bipartite graph which joins an 
unlabelled vertex in X to a labelled vertex in Y, increase the 
current cost by 0. 


Part (c) above will have produced at least one more edge with 
current cost zero. Incorporate all such edges with zero cost in the 
partial graph, delete all labels, and return to Part A, Step 1. 


Summary of algorithm 
START with no matching. 


Assign weights to the vertices and construct the first partial graph. 


Part A: labelling procedure 
Label the vertices to identify an alternating path. 


If none of the vertices on one side of the graph can be labelled, STOP: the 
current assignment is an optimum assignment. 


If breakthrough is achieved, go to Part B. 
If breakthrough is not achieved, go to Part C. 


Part B: matching improvement procedure 
Find an alternating path by tracing back through the labels. 
Form anew matching. 


Return to Part A. 


Part C: modification of partial graph procedure 

Construct a revised cost matrix as follows. 

On the existing cost matrix: 

draw a horizontal line through each labelled vertex in X; 

draw a vertical line through each labelled vertex in Y; 

e find the smallest non-zero entry 6 with ONLY a horizontal line through 
it; 

e decrease all entries with ONLY a horizontal line through them by 4; 
increase the weight on the corresponding vertices in X by 6; 


e increase all entries with ONLY a vertical line through them by 6; 
decrease the weight on the corresponding vertices in Y by 6. 


Construct a revised partial graph. 


Return to Part A. 


Hungarian algorithm for the transportation problem 


START with no flow in a bipartite graph whose vertices are divided into 
two sets: supply vertices Ay, ..., A, and demand vertices by, ..., Dy. 


Construction of the initial partial graph 


STEP 0a Assign weights to the supply and demand vertices as follows. 
To each supply vertex, assign a weight equal to the lowest cost on 
any edge incident with that vertex. 
Decrease the cost on each edge by the weight on the supply vertex 
with which it is incident. 
Repeat this for each demand vertex. 

STEP Ob Construct a partial graph consisting of the vertices of the original 
bipartite graph, together with only those edges whose current cost 
is Zero. 


Part A: labelling procedure 


STEP 1 Label with (*) each supply vertex whose supply has not all been 
allocated to edges. 

If no such vertex exists, STOP: the current solution is a minimum- 
cost solution. 
Otherwise, go to Step 2. 

STEP 2 Choose a newly labelled supply vertex Aj, and label with (A) all 
unlabelled demand vertices which are joined by an edge of the 
partial graph. 

Repeat this for all newly labelled supply vertices, and then go to 
Step 3. 

STEP 3. Choose a newly labelled demand vertex b;, and label with (6;) all 
unlabelled supply vertices which are joined to b; by an edge which 
has been allocated a flow. “STOP: 


Repeat this for all newly labelled demand vertices. 
Repeat Steps 2 and 3 until 


EITHER a demand vertex whose demand is not satisfied is labelled (this is 
called breakthrough), in which case go to Part B, Step 4, 


OR no more labelling is possible, in which case go to Part C, Step 6. 


Part B: flow-augmenting procedure 


STEP 4 Start with the breakthrough vertex, and trace back through the 
labels until a supply vertex which is labelled (+) is reached. A 


flow-augmenting path has been found. The flow-augmenting path plays a 
similar role in this algorithm to that of 
the alternating path in the algorithm 
for the assignment problem. 


Calculate the maximum flow F which can be sent along this path 
as follows. 


First, find the minimum of all flows assigned to EVEN edges of the 

path (that is, the 2nd, 4th, ... edges). This is the maximum possible 

backward flow. 

Next, find the minimum value of: Note that for a flow-augmenting path 
with just one edge, there is no 
backward flow, so we just take the 

e the required demand at the end vertex; minimum of the supply and the 
demand. 


e the available supply at the start vertex of the path; 


e the maximum backward flow. 
This is the required value of F. 
STEP 5 Form anew flow pattern as follows. 


(i) Increase the flows on the ODD edges of the path by this 
value of flow, F. 


Decrease the flows on the EVEN edges by F. 
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(ii) Decrease the available supply at the START vertex by F. 
Decrease the demand at the END vertex by F. 
Remove all the labels and return to Part A, Step 1. 


Part C: modification of partial graph procedure 
STEP 6 Find the lowest current cost on any edge of the original bipartite 
graph which 
e starts ata LABELLED supply vertex; 
e ends at an UNLABELLED demand vertex; 
e has a non-zero current cost. 
Call this lowest cost 6. 
STEP 7 (a) Increase the weight on each labelled supply vertex by 6. 
(b) Decrease the weight on each labelled demand vertex by 0. 


(c) For each edge which joins a labelled supply vertex to an 
unlabelled demand vertex, decrease the current cost by 6. 


(d) For each edge which joins an unlabelled supply vertex to a 
labelled demand vertex, increase the current cost by 0. 


Incorporate all edges which now have a current cost of zero in 
the partial graph, delete all labels, and return to Part A, 
Step 1. 


Summary of algorithm 


START with no flow. 
Construct the initial partial graph. 


Part A: labelling procedure 
Label the vertices to identify a flow-augmenting path. 


If no labelling is possible, STOP: the current solution is a minimum-cost 
solution. 


If breakthrough is achieved, go to Part B. 
If breakthrough is not achieved, go to Part C. 


Part B: flow-augmenting procedure 

Find a flow-augmenting path by tracing back through the labels. 
Augment the flow. 

Return to Part A. 


Part C: modification of the partial graph procedure 

Construct a new revised cost matrix as follows. 

On the existing cost matrix: 

draw a horizontal line through each labelled supply vertex; 

draw a vertical line through each labelled demand vertex: 

e find the smallest non-zero entry 5 with only a horizontal line through 
it; 

e decrease all entries with only a horizontal line through them by 9; 
increase the weight on the corresponding supply vertices by 9; 


increase all entries with only a vertical line through them by 9; 
decrease the weight on the corresponding demand vertices by 0. 


Construct a revised partial graph. 
Return to Part A. 
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Design 3 Design of codes 


Section 1. Error detection and correction 


1 A code is a set of binary words, each containing the 
same number of bits (binary digits 0 and 1). Each binary 
word is a codeword. Such a code is known as a block 
code, because each codeword is a block of binary digits. 
The number n of bits in each codeword of a block code is 
the length of the code. 


2 Most of the codes studied in this unit are constructed 
by adding n — k check bits to each message of length k to 
make a binary word of length n. Such a code consists of 
2* of the 2” possible strings of n bits. 

A code of length n, with k message bits, is an (n, k) code. 
The number k is the dimension of the code. 

The rate of a code is the number of message bits in each 
codeword divided by the length of the code. The rate of 
an (n, k) code is k/ n. 


3 Encoding means identifying each possible binary 
message of (say) k bits with a binary codeword containing 
n bits, where n > k. 

Decoding means estimating the original message from 
the received word. If few errors have occurred in 
transmission and the code is efficient, this can be done 
accurately. 


4 An encoding rule that results in the k bits of the 
message appearing in fixed positions in the 
corresponding codeword is a systematic encoding rule. 
The k message bits need not necessarily form the first k 
bits of each codeword, nor need they appear in k 
consecutive places, but they must appear in the same 
order in each codeword, in some fixed set of k places. 


5 The weight w(x) of a binary word x is the number of 1s 
in x. 
The even-weight code of length 4 is the code 

{0000, 0011, 0101, 0110, 1001, 1010, 1100, 1111}. 


The fourth bit is an overall parity check. The term even- 
weight means that the sum of the bits in each codeword 
is even. 


6 Majority-logic decoding is used when there are 
several ‘votes’ for each bit of information; the value with 
the largest number of votes is chosen. 


7 Repetition codes 
The n-fold repetition code R(n) encodes the messages 0 
and 1 as follows: 


message codeword 
0 00 ... 0 (n bits) 
1 11... 1 (n bits) 


The first of these words is the zero word, denoted by 0. 
The rate of the code R(n) is 1/n. 


8 The Hamming distance d(x, y) between two binary 
words x and y of the same length is the number of places 
in which their bits differ. 

The Hamming distance possesses the following 
properties. 

For any binary words x, y and z of the same length: 

e d(x, y)20; 

e d(x,y)=0 ifand only ifx =y; 

° d(x, y)=dy, x); (symmetry) 

e d(x, z)< d(x, y)+ dy, z). (triangle inequality) 


9 The minimum distance 6 of a code is the smallest 
Hamming distance between two distinct codewords. 

An (n, k) code with minimum distance 6 is sometimes 
referred to as an (n, k, 6) code. 
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Theorem 1.1 
| Let C be a code with minimum distance 6. 
If 6 is odd, then C can detect and correct up to | 
(O-1)/2 errors. 
If 6 is even, then C can detect up to 6/2 errors and 
correct up to (6-2)/2 errors. 


A code with minimum distance 6 can correct | (6- 1)/2] 
errors, where |x | denotes the integer part of x. 


Section 2. Linear codes 


1 A linear code is a code in which the sum (modulo 2) of 
any two codewords is also a codeword — that is, 
whenever x and y are codewords, then so is 
x+y. 


2 The weight w(x) of a binary word x is the number of 1s 
in x. 


Theorem 2.1 
For a linear code, the weight w possesses the 
properties: 
(a) for each pair of codewords x and y, 

d(x, y) = w(x + y); 
(b) the minimum distance 6 of the code is the 
smallest non-zero weight among the codewords. 


Other useful properties are: 
¢ w(x +y)< w(x) + Wy); 
¢ w(x +y)2w(x)- w(y). 


3 Generator sets and matrices 

A generator set for a code is a minimal set of codewords 
with the property that every codeword can be obtained 
by adding codewords in this minimal set. A generator set 
is denoted by angled brackets. 

For a code with dimension k, any generator set consists of 
k codewords, none of which can be formed by adding any 
of the others in the set. 

A code can have several different generator sets. 

Ak xn matrix G is a generator matrix for an (n, k) linear 
code C if the binary words that can be expressed as a 
sum of a subset of the rows of G are exactly the 
codewords of C. 

A code can have several different generator matrices. 

A generator matrix G can be used to specify a simple 
rule for assigning a unique codeword to each possible 
message of k bits. 

Let a message m be a binary word of length k written as a 
row vector. The matrix product mG is a binary word of 
length n, which must be a codeword since mG is a sum of 
rows of G. If the bits mj; , mj, , ..., mj, of m are 1 and all the 
other bits of m are 0, then mG can be obtained by 
adding rows j,,j2,--., jp of G. It can be shown that 2* 
different binary words of length n can be obtained by 
adding rows of G, and hence that no two distinct 
messages m of length k can give rise to the same 
codeword mG. 

When calculating the matrix product mG, addition and 
multiplication modulo 2 are used. The rules are 
summarized as follows. 


+10. x|0 1 
+o 7 0/0 0 
Po | 


Encoding using a k x n generator matrix G is systematic 
whenever k of the n columns of G are distinct columns of 
the k x k identity matrix I, . If column i of Ix forms column 
jj of G, then the ith bit of each message appears as the j; 
th bit of the corresponding codeword. A k x n generator 
matrix is in standard form when the first k columns form 
the k x k identity matrix I,.. 


4 Parity check matrices 

Let C be an (n, k) linear code. 

_ A parity check matrix H is an (n — k) x n matrix with the 
property that the codewords of C are precisely the binary 
words x that satisfy the matrix equation Hx’ = 07. 

A code can have several different parity check matrices. 
A binary word x satisfies Hx’ = 0! if and only if the sums 
(obtained by ordinary addition) of certain of its bits 
(determined by H) are even — that is, they have even 
parity. The n —k linear equations are parity check 
equations of the code. 

Suppose that a codeword x from a code with parity check 
matrix H is transmitted over a communication channel. 
Let e denote the error word that affects x, so that the 
word r =x + e is received. Then 


Hr? =H(x+e)! = H(x! +e7)=Hx! +He! =0! + 
He! 
='He’ . 


In other words, the column vector of length n — k that 
results from premultiplying the transpose of the received 
word r by a parity check matrix H is independent of the 
codeword transmitted. It depends only on the errors that 
have occurred. It is the sum of those columns of H that 
correspond to the corrupted bits. 

For an (n, k) code with parity check matrix H, the error 
syndrome of a received word r is the column vector Hr! 
of length n-k. 


5 Hamming codes 

Consider an (n, k) code that corrects a single error, and 
let m= n-—k. Then a parity check matrix of this code has 
2™ — 1 possible columns. Let H be the matrix formed 
from these columns, written in increasing binary order. 
Then H is a parity check matrix of a Hamming code. 
This code has a simple decoding algorithm. If a 
transmitted codeword is received with its ith bit in error, 
then the syndrome calculated from the received word is 
the ith column of H. However, reading from top to 
bottom, the ith column of H is the binary representation 
of the number i. The bit that is in error can therefore be 
corrected immediately. 

It is possible to construct a Hamming code of length 
' 2™ — 1 for each integer m, where m 2 2. A parity check 
matrix has m rows, and 2” —1columns. The code has 
2™ —1-—m message bits. All Hamming codes have 
minimum distance 3, so each is a (2” -—1,2" —1-m,3) 
code, for some integer m. 


6 Perfect codes 

Suppose that c is a codeword in an (n, k, 5) code, and that 
t is an integer such that t < (6-1)/2 if 6 is odd, and t < (6- 
2)/2 if dis even. Then there are at least 


{eth a) 


n-bit binary words. However, the total number of n-bit 
binary words is 2”, and so 


1+()+(t)e--+()<2" 


This inequality is called the Hamming inequality; the 
upper bound given by this inequality is the Hamming 
bound. 

A perfect code is a linear code for which the Hamming 
bound is attained. 

Every Hamming code is a perfect code. 


Section 3 New codes from old 


1 Equivalent codes 


Two codes are equivalent if the codewords of one can be 
obtained by rearranging the bits of each codeword of the 
other according to some fixed rule. 


Theorem 3.1 
Any (2 -1,2" -—1- m, 3) linear code is equivalent 


to the Hamming code of the same length and 
dimension. 


For this reason, any (2 -—1,2™ -—1- m, 3) linear code is 
known as a Hamming code. 


Theorem 3.2: Every (n, k) code is equivalent to [aeede-with-e- 
systematic encoding tute—thatis] a code whose parity check 


matrix can be written in the standard form H=[AII], where A 
is an (n —k)xk matrix of Os and 1s and I is the (un —k) x (n—k) 


identity matrix. 


If a parity check matrix H of a code C is expressed in 
standard form H = [A|I], then the matrix G = [I| A'] is 
a generator for C. [And conversely, if the generator 
matrix G = [I| A’], then the matrix H=[AlI] isa 
parity check matrix for C.] 


A cyclic code is a code of length n with the property that, 
whenever x =X, X ... X, is a codeword, then so also is 
X' =X)... X_ X). For each positive integer m, there is a cyclic 
Hamming code of length 2” - 1. In other words, each 
Hamming code is equivalent to a cyclic Hamming code. 


3 Dual codes 

Let H be a parity check matrix of a linear code C. The 
sum of each possible subset of rows of H is a unique 
binary word. Thus exactly 2"~* different words of length 
n can be obtained from the rows of H. These words are 
therefore the codewords of an (n, n —k) code, called the 
dual code C” of the original code C. 

A simplex code is the dual of a Hamming code: the dual 
of a (2™ -—1, 2" -—1- m, 3)Hamming code is a 
(27 —1,m,2™~ 1) simplex code. 

The simplex code of length 2"-—1 contains the zero 
word and 2” — 1 words of weight 2”~ ' . If the codewords 
are considered as vertices of the n-cube, they form a 
regular n-simplex — that is, a simplex with edges of the 
same length. 

Let H be a parity check matrix of an (n, k) linear code C. 
The dual of C is the (n, n —k) linear code C” with 
generator matrix H. 


Theorem 3.3: duality theorem 
Let C be a linear code. Then (C” )" = C. 


Corollary 


Let C be a linear code. Then a parity check matrix of 
the dual code C” is a generator matrix for C. 


Theorem 3.4 
The dual of a cyclic code is cyclic. 
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4 Extended codes 

A given (n,k, 6) code can be extended by adding an 
overall parity check. An extra bsit is added to each 
codeword in such a way that the resulting binary word of 
length n + 1 has even weight, so a 0 is added to each 
codeword of even weight, and a 1 to each codeword of 
odd weight. The minimum distance of the extended code 
is 6+ 1if 6 isodd, and dif 6 is even. 

A parity check matrix for the code obtained by extending 
an (n, k) code with parity check matrix H is 


Be Cat Bite 
De oe een 
0: H 

0! 


(Here we have added the overall parity check as the first 
bit of each codeword.) 
5 The [ala + b] construction 


Theorem 3.5 


Let A be an (n,k4, 64 ) code and B be an (n,k g, 63 ) 
code, and let 6 be the smaller of 26, and 6, . Then 


the code C whose codewords are all the binary words 
of the form [ala+b], where a is a codeword of code A 
and b is a codeword of code B, is a (2n,k4 +kpg, 6) 
code. 


If G4 is a generator matrix for code A and Gg is a 
generator matrix for code B, then a generator matrix for 
code C is 


If Hy, is a parity check matrix for code A and Hz is a 
parity check matrix for code B, then a parity check matrix 
for code Cis 


Section 4 First-order Reed-Muller codes 


1 Three special codes 

R (1) is the code {00, 01, 10, 11}. 

R (2) is the even-weight code of length 4. 

R (3) is the (8, 4, 4) extended Hamming code. 


2 Reed-Muller codes 
A Boolean variable is a variable that can take only two 
values, usually 0 and 1. A Boolean function is a function 
of Boolean variables, which can itself take only the same 
two values. 
The set of binary words that correspond to all the 
Boolean functions of degree at most 1 in m variables, v 
(051.0 Um , is called the first-order Reed-Muller code o 
length 2™,; itis denoted by R(m). 
The dimension of ® (m) is m+1. 
The code & (m) has: 

one codeword of weight 0, 

2*1 — 2 codewords of weight 2", 

one codeword of weight 2” . 
The minimum distance of R(m)is2™'. 
The code  (m) isa (2”™,m +1,2™~' ) code. 


Theorem 4.1 
The code R(m + 1) consists of all binary words of the 
form [ala+b], where a is a codeword in R(m) and b is. 
a codeword in the repetition code R ( 2” ). 


A generator matrix for ® (3) is 


1: dy: vba dcak. « tile I 
ey a ea. 
aoe flees Site! ote rie aaa fe 
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Let G, be a generator matrix for R(m). Then a generator 
matrix for R(m +1) is 


First-order Reed-Muller codes can be decoded using a 
form of majority logic. 

Suppose a codeword is transmitted that corresponds to 
the message a, encoded using the generator matrix G. 

Let r be the word received. If no error has occurred in 
transmission, then r = aG. This implies relationships 
between the bits of r and a which give several ‘votes’ for 
each bit of a. 


3 Golay codes 
Let A denote the first-order Reed-Muller code & (3) of 
length 8 with generator matrix 


Podangh. Gt: sb 139 
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If the order of the first seven bits of each codeword of 
code A is reversed, the codewords of an equivalent first- 
order Reed-Muller code of length 8 are obtained. 
Denote this second code by B. 
The extended Golay code G(24) is the code whose 
codewords are all the binary words that can be written in 
the form 

[aj +bla, +bla, +a, +b], 
where a, and a, are codewords in A, and b is a codeword 
in B. 
The code G(24) is a (24, 12, 8) code. 
Let G, and Gg be generator matrices of codes A and B, 
respectively. Then a generator matrix for the code G(24) 
is the 12 x 24 matrix 


Ga oly rey 
G, ! G, ! B 


The code G(24) is self-dual. 

The Golay code G(23) is a (23, 12, 7) code obtained from 
G(24) by deleting the last bit of each codeword. 

The perfect codes are: the repetition codes R(n) for all 
odd values of n, the Hamming codes, and the two Golay 
codes. There are no others. 


4 Mariner 9 code 
The Mariner 9 code is the code & (5). 


Graphs 4 Graphs and computing 


Section1 Efficiency of algorithms 


1 The time complexity function T(n) for an algorithm 
gives the maximum time taken to process any input of size 
n, at 1 time unit per comparison made. 


2 The space complexity function S(n) for an algorithm 
gives the maximum space required in a computer’s 
memory for storing the data used by the algorithm. 


3 Let T(n) be a time complexity function, and let g(n) be a 
function for which there exists both a positive constant c 
and a (large enough) number N such that 

T(n)<c.gn), foralln2N. 
Then g(n) (asymptotically) dominates T(n), and T(n) is 
(asymptotically) dominated by g(n). 


The set of all functions that g(n) dominates is denoted by 
O(g(n)). 


If g(n) dominates T(n), then T(n) is O(g(1)), meaning T(r) is 
in the set O(g(n)). 


If a function g(n) dominates a time complexity function 
T(n) and T(n) also dominates g(n), then T(1) has order 
O(g(n)), and T(n) and g(r) have the same order of 
magnitude. 


4 The order hierarchy of common big-oh sets is 
O(1) < Oflogyn) < O(n) c O(nlogyn) € O(n’) c ... € 
id on’) c ...€ O (2") c O(n!). 
slow 
The time complexity function T(n) of an algorithm is 
placed in the set O(g(n)) in the hierarchy if T(n) has order 


O(g(n)). The further to the left in the hierarchy T(n) is 
placed, the faster the algorithm (for large enough n). 


To find the order of a given time complexity function T(1) 
consisting of the sum of one or more terms: 

e take the base of each logarithm to be 2; 

e take each coefficient to be 1; 

e take each constant term to be 1. 

Then the term g(n) whose set lies furthest to the right in 
the order hierarchy is the dominant term, and T(n) has 
order O(g(1)). 


Section 2. Stacks and lists 


1 The stack data type consists of data stored in an 
ordered set (stack) of adjacent storage cells together with 
the basic operations: 

TOP(s) = top item in stack s; 

DEPTH(s ) = number of items in stack s; 

PUSH(item, s) = stack created by pushing (adding) item 
onto top of stack s; 

POP(s) = stack created by popping (removing) top item 
from stack s; 

ISEMPTYSTACK? = TRUE if s is the empty stack, and FALSE 
otherwise. 


2 The list data type is a data store k in which each item 
keeps an address that points to the next item, together 
with the basic operations: 

FIRST(k) = first item in list k; 

ITEM(i, k) = ith item in list k; 

LENGTH(k) = number of items in list k; 

INSERT(item, i, k) = insert item after (i— 1)th item in list k so 
that it becomes the ith item of the new list formed; 
DELETE(i, k) = remove ith item from list k. 


3 A depth-first search of a rooted path graph starts at the 
root vertex and moves from level to level, comparing the 
item at each vertex with the given item. 


Theorem 2.1 


The time complexity function for a depth-first search of 
a rooted path graph is of order O(n). 


Algorithm: SEARCH 1 for a given item in a list k of 

length n 

STEP 1 Set the value of a variable i = 0. 

STEP 2 Increase the value of i by 1 and compare ITEM(/, k) 
with the given item. If they are the same, or if they 
are not the same but i = n, sToP. Otherwise repeat 
Step 2. 


Algorithm: SEARCH 2 for a given item in a list k 

STEP 1 If LENGTH(k) = 0, the list is empty, so STOP. 

STEP 2 Compare the given item with ITEM(1, k). If they are 
the same, STOP. 

STEP 3 Use SEARCH 2 on the list DELETE(1, k). 


4 Arecurrence system is a system of equations of the 
form 

T(k) = a(n) 

T(n) = b(n) + T(n-1) (n >k). 
Solution by iteration is the derivation of a formula for 
T(n) using a recurrence system. 


5 


Theorem 2.2 


The time complexity function for a search for all 
repeated items in a rooted path graph is of order O(n’). 


6 The bubble sort algorithm takes the form of a number 
of ‘passes’ over all the items in a list. For the first pass, 
start at the first item in the list, at the root vertex, and 
compare it with the item at the second vertex. If the first 
item is ‘larger’, exchange the two vertices; if the first item 
is not larger, do nothing. Work along the graph, 
comparing the ith vertex with the (i+ 1)th vertex, for each 
i, and exchanging them if necessary. The number of such 
passes made is n - 1, each pass using one fewer 
comparison than the previous one. 


Theorem 2.3 


The time complexity function for a bubble sort is of 
order O(n’). 


7 Algorithm to merge two sorted lists 

STEP 1 Start a new list, with no items in it. 

STEP 2 Compare the first items in both sorted lists, delete 
the smaller of these from its list, and insert it at the 
end of the new list. Repeat until one list is empty. 

STEP 3 Insert what remains of the non-empty list at the 
end of the new list, and STOP. 


Theorem 2.4 


The time complexity function for merging two sorted 
lists of lengths m and n, where m < n, is of order O(n). 


8 The merge sort algorithm first bisects a list into two 
lists of equal (or nearly equal) length. It then bisects each 
such half-length list into two equal (or nearly equal) 
quarter-length lists. The bisection process continues until 
a number of pairs of short lists are obtained, each with one 
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or two items. These short lists are easily sorted, as there is 
nothing to do for a one-item list, and a single comparison 
to make for a two-item list. The pairs of short lists are now 
merged to give pairs of lists, each with three or four items. 
These pairs are merged, and so on, until the whole list is 
sorted. 


Theorem 2.5 


The time complexity function for a merge sort on a list 
_ of length n is of order O(nlog, n). 


9 The divide-and-conquer method of problem solving 
breaks down a problem into subproblems, solves these 
subproblems, and then combines the solutions to the 
subproblems into a solution for the original problem. The 
method is recursive; the subproblems themselves can be 
solved using the divide-and-conquer method. 


Section 3 Binary trees 


1 For each vertex v of a binary tree, an adjacent vertex 
below and to the left can be regarded as the root of 
another binary tree, the left subtree of v, and similarly an 
adjacent vertex below and to the right is the root of the 
right subtree of v. 


2 Algorithm: GROW-TREE for constructing a binary tree 

from a list 

STEP 1 If the list is empty, STOP. 

STEP 2 Find the centre or bicentre of the graph of the list. 
If it has a centre, choose this as the root of the 
binary tree; if it has a bicentre, choose the 
rightmost vertex of the two as the root. 


STEP 3 The root divides the list into a left list and a right 
list (one or both of which may be empty). Apply 
GROW-TREE to both of these to obtain the left and 
right subtrees. 


The use of the GROW-TREE algorithm on lists of lengths 
1,..., 7 yields the following catalogue of binary trees. 


a 
1-tree 2-tree 3-tree 
4-tree 5-tree 
6-tree 7-tree 


3 The height of a binary tree is the number of vertices in 
a path of maximum length, starting from the root. 


A binary tree is balanced if, at each vertex v, the heights of 
the left subtree of v and the right subtree of v do not differ 
by more than 1. 
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a 
4 The binary tree data type is a data store T, 
corresponding to a binary tree, in which each item keeps 
two addresses that point to its two subtrees, together with 
the basic operations: 
ROOT(T) = item at root vertex; 
LEFT(T) = left subtree of T; 
RIGHT(T) = right subtree of T; 
MAKE-TREE(T, , item, T, ) = tree constructed with the item 
as its root, T; as its left subtree and T, as its right subtree; 
ISEMPTYTREE? = TRUE if T is the empty tree, and FALSE 
otherwise. 


5 Given a set of items that can be ordered, a binary 
search tree is a binary tree in which each vertex v 
represents one of these items in such a way that 

e v is larger than each item below and to the left of v; 

e v is smaller than each item below and to the right of v. 


Algorithm: SEARCH a binary search tree 

STEP 1 Compare the search item with the ROOT. If they are 
the same, STOP. 

STEP 2 If the search item is smaller than the ROOT, then 
apply SEARCH to the LEFT tree, if it exists. Otherwise 
Srur. 


If the search item is larger than the ROOT, then 
apply SEARCH to the RIGHT tree, if it exists. 
Otherwise STOP. 


Theorem 3.1 | 


The time complexity function for a search on a binary | 
search tree of height / is of order O(/). 


The time complexity function for a search on a 
balanced binary search tree with n vertices is of order 
O(log, 7). 


6 The idea of depth-first search of a graph is to penetrate 
as deeply as possible into the graph before fanning out to 
the other vertices. 


Algorithm: DEPTH-FIRST SEARCH on a binary tree 
STEP 1 Start at the ROOT. Add it as the rightmost vertex of 
the list. If there are no subtrees, STOP. 


STEP 2 If there is a LEFT subtree, then apply DEPTH-FIRST 
SEARCH to the LEFT subtree. 


If there is a RIGHT subtree, then apply DEPTH-FIRST 
SEARCH to the RIGHT subtree. 


Algorithm for manipulating the stack in a depth-first 
search 
STEP 1 PUSH the root vertex onto the empty stack. 


STEP 2 If the top vertex of the stack is adjacent to a new 
vertex, PUSH this vertex onto the stack. Otherwise, 
POP the top vertex off the stack. 


STEP 3 If the stack is empty, STOP. Otherwise, return to 
Step 2. 


This algorithm can be applied to any rooted tree, not just 
to a binary tree. The algorithm can be used to search any 
connected graph: the only change is that the ‘root vertex’ 
in Step 1 becomes the ‘start vertex’ — any vertex can be 
chosen as the start vertex. 


The tree consisting of the vertices visited and the edges 
traversed in a depth-first search on a connected graph is 
called a depth-first search spanning tree. 


7 A queueis a data store similar to a stack; the difference 
is that a queue is FIFO (First In, First Out) whereas a stack 
is LIFO (Last In, First Out). 


8 The idea of breadth-first search of a graph is to fan out 
to as many vertices as possible before penetrating deep 
into the graph. This means that all the vertices adjacent to 
the current vertex must be visited before going on to 
another vertex. 


Algorithm for manipulating the queue in a breadth-first 
search 
STEP 1 INSERT the root vertex into the empty queue. 


STEP 2 If the first vertex of the queue is adjacent to a new 
vertex, INSERT this vertex at the end of the queue. 
Otherwise, DELETE the first vertex of the queue. 

STEP 3 If the queue is empty, sTOP. Otherwise, return to 
Step 2. 


The tree consisting of the vertices visited and edges 
traversed in a breadth-first search on a connected graph is 
called a breadth-first search spanning tree. 


S 


Theorem 3.2 
The time complexity functions for depth-first search 


and breadth-first search of a connected graph with 
m edges are both of order O(m). 


Section 4 Quad trees 


1 A k-screen consists of 2"! x 2* pixels (small squares). 
p q 


A quad tree that represents a k-screen is a rooted tree with 
at most k levels. Each end-vertex has 1 or 0 stored at it, 
denoting that the pixel (or quadrant of pixels) it represents 
is black or white, respectively. Each other vertex has a set 
of four vertices below it, ordered from left to right, this 
ordering being denoted by the labels A, B, C, D. At each 
level, each set of four vertices corresponds to four 
quadrants of the screen — the A-vertex to the top-right 
quadrant, the B-vertex to the top-left quadrant, and so on, 
in an anticlockwise direction. 


2 To rotate an image on a k-screen anticlockwise through 

a right angle, rotate it about its centre. This rotation is 

represented by a cycle of the quadrants 
A>B-C-»D-A. 

Other rotations can be represented similarly. 


To reflect an image in a horizontal line through the middle 
of a k-screen, interchange the quadrants A and D and the 
quadrants B and C, according to the scheme 

AD and Bec. 
Other reflections can be represented similarly. 


To magnify an image in one of the principal quadrants of 
a k-screen by a factor 2, interpret the subtree for that 
quadrant as the quad tree for the whole screen. 


Magnification by other powers of 2 is also possible by 
moving further down the levels of the quad tree and 
detaching the subtrees. 


To reduce an image on a k-screen by factors of 2, move all 
vertices down one or more levels. 


3 ‘The pixel immediately above another pixel is called the 
north neighbour of the pixel. 


For (blocks of) pixels, the B-quadrant and C-quadrant are 
opposites, and the A-quadrant and D-quadrant are 
opposites, as far as being north neighbours is concerned. 


Algorithm for finding the north neighbour of a vertex of 

a quad tree 

STEP 1 Set up an empty ‘path’ stack to keep track of the 
vertices that are passed. Make the start vertex the 
‘current’ vertex. 


STEP 2 If the ‘current’ vertex is the root vertex, and a 
lower (C or D) vertex is not on the Top of the ‘path’ 
stack, then no north neighbour exists, so STOP. 


Otherwise, PUSH the ‘current’ vertex label onto the 
‘path’ stack. Move up to the parent vertex and 
make it the ‘current’ vertex. If the ToP of the ‘path’ 
stack is not a lower (C or D) vertex, repeat Step 2. 


STEP 3 POP a vertex label from the ‘path’ stack. Move 
down to the vertex with the opposite label to the 
one just taken off the stack, and make it the 
‘current’ vertex. 

Repeat Step 3 until the ‘path’ stack is empty. The 
‘current’ vertex is then the north neighbour of the 
start vertex. 

The north neighbour algorithm can be adapted to yield 

(the colours of) the other seven neighbouring pixels — the 

south, east, west, north-east, north-west, south-east and 

south-west neighbour pixels. 

The algorithms for finding one, or all, of the eight 

neighbours of a given pixel are all of order O(k). 


Time complexity functions for computer algorithms on n items of data 
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Section5 Branch-and-bound methods 

1 A state space tree for a problem is a rooted tree in 

which the root vertex corresponds to the problem as a 

whole and the other vertices correspond to subproblems. 

The subproblems at the end-vertices are ones that can be 

solved immediately. Each parent vertex is: 

either an AND vertex, in which case its (sub)problem can 
be solved immediately once all the subproblems at 
its children vertices have been solved; 

or an OR vertex, in which case its (sub)problem can be 
solved immediately once any one of the 
subproblems at its children vertices has been 
solved. 

An AND tree is a state space tree all of whose parent 

vertices are AND vertices. 

An OR tree is a state space tree all of whose parent vertices 

are OR vertices. 

And AND/OR tree is a state space tree that has a mixture of 

parent vertex types. 
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2 Branch-and-bound method 

Given a problem that requires the optimization of some 
objective, determine an appropriate bound for the 
objective. Break the problem down into a set of 
subproblems and determine the bound for each. 


For each subproblem that is not immediately solved or 
cannot immediately be shown to have a non-optimal or no 
solution, break it down further into subproblems and 
determine the corresponding bounds. Repeat this process 
until each subproblem is solved or has been shown to have 
a non-optimal or no solution. 

An optimal solution to the original problem is given by 
any subproblem with a ‘best’ bound. 


3. The divide-and-conquer method corresponds to an AND 
tree, since every subproblem has to be solved in order to 
solve the initial problem. The branch-and-bound tree is an 
OR tree, so that the solution to each subproblem is a 
solution to the original problem. 


4 Knapsack problem 

A hiker is planning a journey, but has a knapsack that can 
accommodate only a certain total weight. There are a 
number of items that the hiker wishes to take, each of 
which has a particular value. Which items should be 
packed so that the total value of the packed items is a 
maximum, subject to the weight restriction? 


For a problem with n items x, ..., X,, a solution vector is a 
vector of the form (x1, ..., X;,), 


x; =1 if item? is packed; 


x; =0 if item is not packed. 


where 


A feasible solution is one which satisfies the weight 
constraint. 


An algorithm for the knapsack problem 


START with the zero vector (0, 0, ..., 0); this is feasible 
with value 0. 
STORE (0, 0, ..., 0) and value 0. 
GENERAL  e_ Branch from the first solution from 
STEP which branching is possible. 
e Calculate the total weight of each new 
solution. 
e Calculate the total value of each new 
feasible solution. 
If there is a new feasible solution with 
value greater than the value stored, 
store the new solution vector and its 
value instead. 
e Mark a vertex with a U1 if it corresponds 
to: ; 
O avector which equals or exceeds 
the weight restriction; 
Oa vector which ends in 1. 


REPEAT the GENERAL STEP until no more branching 
is possible. 
STOP The stored solution vector and its value gives 


an optimum solution. 
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5 An algorithm for the travelling salesman problem 


START 


GENERAL 
STEP 


REPEAT 


SLOF 


with a given n x n table of distances, 
corresponding to a complete weighted graph 
with n vertices. 

Carry out the initial row and column 
reduction, and calculate the initial lower bound 
for the total length of the route. 


e Consider all the edges with zero 
weight, and choose an allowable 
edge e whose exclusion leads to the 
greatest increase in the lower bound; 
if there are several such edges, choose 
the first. 

(Avoid cycles with fewer than n edges.) 

e Consider the consequences of including 
e and excluding e. Use row and column 
reduction to determine these 
consequences, in terms of increases in 
the lower bound. Choose the option 
which gives the smaller lower bound; if 
the lower bounds are equal, include the 
edge e. 

STORE the current list of included 
edges, and the current lower bound. 

e Continue from the current position 
unless the chosen option has a lower 
bound greater than a previously 
eliminated option, in which case 
backtrack to the earlier position. 


the GENERAL STEP until a cycle with 
edges has been created. 


The stored list of edges gives an optimum 
solution. 


Networks 4 Physical networks 


Section 1 Modelling physical networks 


1 Through and across variables 

A physical network can be considered to be a number of 
individual components connected together. 

In this unit it is assumed that the behaviour of a 
component does not change if it is moved from one system 
to another, and that the behaviour is time-invariant and 
linear. 


The point at which one component is connected to another 
is called a terminal. 


An across variable represents the difference in a quantity 
between one terminal of a component and another. A 
through variable represents the quantity flowing through 
a component. 


2 Equations relating through and across variables 
Kirchhoff’s current law 

The current flow into any vertex of an electrical network is 
equal to the current flow out of it. | 


_ Kirchhoff’s voltage law 

The algebraic sum of the potential differences across all 
the components around any circuit (cycle) in an electrical 
network is zero. 


Analogues of these laws can be given for many types of 
network containing through and across variables. The 
analogue of Kirchhoff’s current law for a network carrying 
a flow is the flow conservation law: what flows into a vertex 
of a network (other than a source or sink) must flow out of 
it. The analogue of Kirchhoff’s voltage law is the potential 
difference conservation law: the sum of the potential 
differences around any circuit of a network is zero. 


3 Graphical representation of components 

A 2-terminal component can be represented by two 
vertices joined by a single edge. 

Associated with the edge are two edge variables — one 
through variable and one across variable. 

A reference direction for the edge variables is obtained by 
assigning an arbitrary direction to each edge, indicated by 
drawing an arrow alongside the edge. Currents which 
actually flow through the component corresponding to the 
edge in the direction of the arrow are regarded as positive, 
and currents in the opposite direction are regarded as 
negative. Similarly, voltage differences between the 
terminals corresponding to two vertices are always taken 
as the voltage at the tail of the arrow minus the voltage at 
the head. Thus if the voltage at the tail is greater than the 
voltage at the head, the voltage difference is positive: if the 
voltage at the tail is less, then the voltage difference is 
negative. The reference direction is purely arbitrary and 
does not necessarily represent the actual direction of the 
current flow or voltage difference. 


Any system consisting entirely of 2-terminal components 
can be represented by a graph whose edges correspond to 
the components and whose vertices correspond to the 
terminals. Such a graph, with a reference direction 
associated with each edge, is called an oriented graph. 


4 Two-terminal components 
An equation relating the two edge variables for a 
2-terminal component is a component equation. 


A resistor is modelled by a 2-terminal component whose 

edge variables v and i satisfy a relationship of the form 
v=RtL 

where R is a constant, usually called the resistance. 


This relationship is known as Ohm’s law. It is sometimes 
written in the form i= Gv, where G (= 1/R) is the 
conductance of the resistor. 


A capacitor is modelled by a 2-terminal component whose 
edge variables v and i satisfy a relationship of the form 
_C dv 
crea 
where C is a constant, usually called the capacitance. 


An inductor is modelled by a 2-terminal component 
whose edge variables v and i satisfy a relationship of the 
form - 
i 
ve=L dt , 
where L is a constant, usually called the inductance. 


Components such as batteries and electric generators are 
modelled by the concept of an ideal independent source. 
An ideal independent voltage source maintains a prescribed 
voltage across its terminals, irrespective of the size of the 
current flowing between them. An ideal independent 
current source maintains a prescribed current, irrespective of 
the value of the voltage across its terminals. 


Some important types of component, together with their 
graphical representations and their corresponding 
component equations are listed on page 42. 


5 Multi-terminal components 


A multi-terminal component is a component which has 
more than two terminals, and which cannot therefore be 
represented by a single edge joining two vertices. Two 
important types of electrical multi-terminal component are 
transformers and transistors. 


Graphical representations and component equations of 
some multi-terminal components are given on page 43. 


A 3-terminal component is represented by two edges and 
has two component equations. 

Similarly, any labelled tree with n vertices is a suitable 
representation of an n-terminal component. It follows that 
any n-terminal component can be regarded as n—1 2-terminal 
components, one corresponding to each edge of the tree. 


6 Kirchhoff’s laws (alternative forms) 


Kirchhoff’s laws can be formulated in graph-theoretical 
terms as follows. 


Vertex law 


The algebraic sum of the through variables associated with 
the edges at each vertex of an oriented graph is zero. 


Here ‘the algebraic sum’ means that the orientation of the 
edges must be taken into account, so if edges oriented into 
a vertex are counted as positive, then edges oriented out of a 
vertex are counted as negative. 

Cycle law 

The algebraic sum of the across variables associated with 
the edges in any cycle of an oriented graph is zero. 

Here ‘the algebraic sum’ means that a particular direction 
around a cycle (say clockwise) is chosen, and edges 
oriented in the same direction are considered as positive and 
edges oriented in the opposite direction as negative. 


7 Dual electrical networks 


The idea of duality can be applied to (planar) electrical 
networks as follows. 

Given an electrical network N, draw a plane drawing of a 
corresponding oriented graph, form the dual of this planar 
graph as described in Graphs 3, and then construct the 
corresponding network N*, called the dual network. 
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The rule for associating an arrow with each edge of the 
dual graph is illustrated below. 


edge of G* 
edge of G 


go clockwise from edge of G 
to edge of G* in same direction 


Voltage in the original network N corresponds to current in 
the dual network N* and vice versa; the roles of v and iin 
the component equations are interchanged. The following 
table shows correspondences between the components in 
the two networks N and N*. 


network N dual network N* 
voltage current 
current voltage 


resistor (resistance 1/R) 
inductor (inductance C) 
capacitor (capacitance L) 


resistor (resistance R) 
capacitor (capacitance C) 
inductor (inductance L) 


Section 2 Electrical networks: 
matrix equations 


1 Linearly independent equations 

A number of equations are linearly dependent if at least 
one of them ‘depends’ on the others, in the sense that it 
can be obtained by adding or subtracting multiples of the 
others. If none of the equations depends on the others in 
this way, then the equations are linearly independent. 


2 Fundamental cycles 

A number of cycles in a graph are linearly independent if 
the corresponding voltage equations are linearly 
independent. 


Let G be a connected graph. A spanning tree in G is a 
subgraph of G which includes all the vertices of G and is 
also a tree. The edges of the tree are called branches. If all 
the edges of a spanning tree are removed from G, the 
remaining subgraph is called a co-tree and its edges are 
called chords. 


Let G be a connected graph, and let T be a spanning tree of 
G. The set of fundamental cycles associated with T 
consists of the cycles of G each of which is obtained by 
adding a single chord to T. The corresponding voltage law 
equations are called the fundamental cycle equations. 


If G has n vertices and m edges, then any spanning tree 
gives rise to m—n-+1 chords, and so there are m-—n+ 1 
cycles in a fundamental set. This number is the largest 
number of linearly independent cycles in G, and is called 
the cycle rank of G. 


In any electrical network problem, all the information 

contained in Kirchhoff’s voltage law can be obtained by: 

(a) choosing a spanning tree; 

(b) finding the fundamental cycles associated with it; 

(c) finding the corresponding fundamental cycle 
equations. 


The fundamental cycle equations can be represented in 
matrix form. The rows of the matrix correspond to the 
chords associated with a spanning tree, and the columns 
correspond to all the edges of the graph G. (It is usual to 
write the columns corresponding to the branches of the 
spanning tree first, followed by those for the chords.) If G 
has n vertices and m edges, an (m —n + 1) Xm matrix is 
obtained. This matrix is called the fundamental cycle 
matrix, and is denoted by C,. 
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To construct this matrix, add each chord to the spanning 
tree in turn, and look at the edges in the corresponding 
fundamental cycle. Trace around this cycle in the direction 
of the chord, then fill in the row corresponding to the 
given chord, writing: 
1 ineach column corresponding to an edge of the 
cycle oriented in the same direction as the chord; 
-1 ineach column corresponding to an edge 
oriented in the opposite direction; 
0 ineach column corresponding to an edge which 
is not in the cycle. 


The fundamental cycle equations are written in matrix 
form as 

C; v=0, 
where v is the column vector of edge voltages, written in 
the same order as the columns of the matrix, and 0 is the 
appropriate zero vector. 


3 Fundamental cutsets 

Let G be a connected graph, and let T be a spanning tree of 
G. The set of fundamental cutsets associated with T 
consists of the cutsets of G obtained by removing a branch 
of T, thus separating the tree into two parts X and Y, and 
listing the edges of G joining a vertex in X and a vertex in 
Y. The corresponding current law equations are called the 
fundamental cutset equations. 


A number of cutsets in a graph G are linearly 
independent if the corresponding current law equations 
are linearly independent. Thus the cutsets in a set of 
fundamental cutsets are linearly independent, since each 
cutset contains one edge (the branch of the spanning tree) 
not contained in any of the other cutsets. If G has n 
vertices, then any spanning tree has n —- 1 branches, and so 
there are n — 1 cutsets in a fundamental set. This number is 
the largest number of independent cutsets in G, and is 
called the cutset rank of G. 


In any electrical network problem, all the information 

contained in Kirchhoff’s current law can be obtained by: 

(a) choosing a spanning tree; 

(b) finding the fundamental cutsets associated with it; 

(c) finding the corresponding fundamental cutset 
equations. 


The fundamental cutset equations can be represented in 
matrix form. The rows of the matrix correspond to the 
branches of a spanning tree, and the columns correspond 
to all the edges of the graph G. (It is usual to write the 
columns corresponding to the branches first, followed by 
those for the chords.) If G has n vertices and m edges, an 
(n — 1) X m matrix is obtained. This matrix is called the 
fundamental cutset matrix, and is denoted by Dy. 
To construct this matrix, take each branch of the spanning 
tree in turn, and look at the edges in the corresponding 
fundamental cutset. Then fill in the row corresponding to 
the given branch, writing: 
1 ineach column corresponding to an edge of the 
cutset oriented from X to Y; 
-1_ ineach column corresponding to an edge of the 
cutset oriented from Y to X; 
0  ineach column corresponding to an edge not in 
the cutset. 


The fundamental cutset equations are written in matrix 
form as 

D; i= 0, 
where i is the column vector of edge currents, written in 
the same order as the columns of the matrix, and 0 is the 
appropriate zero vector. 


4 Obtaining the fundamental cycle and cutset matrices 
Let G be an oriented graph with n vertices and m edges. 
The incidence matrix B(G) is the n x m matrix in which the 
entry in row 1 and column is 
1 if edge jis incident with, and oriented away 
from, vertex 1; 
-1 if edge j is incident with, and oriented towards, 
vertex 1; 
0 if edge] is not incident with vertex 1. 


Any single row of B(G) can be omitted without loss of 
information. The vertex corresponding to the row omitted 
is the reference vertex and the resulting matrix is the 
reduced incidence matrix, denoted by Bp. 


The matrix By is partitioned into two parts, one part 
corresponding to the branches, and one part 
corresponding to the chords. This partitioned matrix can 
be written as B, = [B, | B,J], where B, is the tree part and B, 
is the co-tree part. 


Theorem 2.1 


The fundamental cycle matrix C; and the fundamental 
cutset matrix D; can be expressed in terms of B, and B, 
as follows: 


fundamental cycle matrix C; = |-(B;"B. : | i et 


fundamental cutset matrix D, = [I ae | B,'B . ; 


5 Tellegen’s theorem 


Theorem 2.2: Tellegen’s theorem 
Suppose that there are two electrical networks which 
can be represented by the same oriented graph with n 
vertices and m edges. Then 


™m 
SY vel = Vi + Vol Sie Diskin =, 

k=1 
where 1; is the voltage associated with the kth edge for 
the first network, and i, is the current associated with 
the kth edge for the second network. 


Section 3 Electrical networks: 
solving the network equations 


1 Combining the matrix equations 

There are three matrix equations, corresponding to the 
component equations, the fundamental cycles, and the 
fundamental cutsets. | 

When the voltages and currents in each equation are 
written in the same order, the three matrix equations can 
be combined into one as follows: 


(component equation) 


(fundamental cycle equation) 


(fundamental cutset equation) 


which can be written as one equation 


H_ isa2mxX2m matrix containing component, 
cycle and cutset information; 

x is acolumn vector containing all the voltages 
and currents; 

y isacolumn vector consisting mainly of zeros, 
but also including terms such as I and —V 
(source functions). 


2 Solving the matrix equation Hx = y } 
Gaussian elimination .is a systematic way of solving a 
linearly independent system of n simultaneous linear 
equations in n variables by eliminating the variables one at 
a time. This is done by performing certain allowable 
operations on the rows of the matrix of coefficients H so as 
to reduce all the elements below the main diagonal (the 
diagonal from top left to bottom right) to zero. The 
equations Hx = y can then be successively solved, starting 
with the last equation, until all the required variables are 
obtained. 


The operations allowed on the rows of the matrix are: 

e interchange two rows; 

¢ multiply any row by a non-zero number; 

¢ add any row to any other row. 

If these row operations are also carried out on the column 
vector y, then each of them corresponds to an operation on 
the corresponding equations which leaves the solution of 
the equations unchanged. 


3 State equations 
For networks involving capacitors and inductors, the 
component equations involve not only v and 1, but also 
their derivatives. The matrix equation can be adapted to 
take the derivatives into account. 
The matrix equation can be written in the form 

Hx = y + Kx, 
where 

x is a column vector involving all the voltages 

and currents; 


x is acolumn vector involving their derivatives; 


y is the ‘source function vector’ introduced 
earlier; 
H_ is the matrix involving the components and the 
fundamental cycles and cutsets; 
“Kis a matrix isolating the derivatives of the edge 
variables. 


The variables whose derivatives occur in the component 
equations appear Iast of all in the column vectors x and x. 
These variables are called state variables. The values of 
the state variables represent the state of the system. 


The state equations are first-order differential equations 
that relate the derivatives of the state variables to the state 
variables and the independent sources. 
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Graphical representation of 2-terminal components 
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component 
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p =Kq 


(The relationship f= mie for a spring is often written in the form f= kx, known as 


Hooke’s law, where x is the extension of the spring due to an applied force f.) 
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Graphical representation of 3-terminal components 
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Design 4 Block designs 


Section 1. Blocking in designed 
experiments 


1 Varieties, plots and blocks 

In an experiment, the varieties are the items being 
compared. 

A plot is a unit of experimental material on which a single 
variety is tested. 


A block is a group of plots which are similar in the sense 
that, if the same variety is tested on each, then similar 
results can be expected. 


2 Incomplete block designs 

An incomplete block design consists of: 

e aset of varieties; 

e aset of blocks, each containing the same number of 
varieties; 


e a list showing the varieties allocated to the plots in each 
block. 


No block contains every variety, and no block contains 
any variety more than once. 


Notation 

v denotes the number of varieties; 

b denotes the number of blocks; 

k denotes the block size — the number of plots in each 
block. 


A block design is equi-replicate if each variety occurs the 
same number of times; this number is the replication, and 
is denoted by r. 


In this unit, all incomplete block designs are assumed to be 
equi-replicate, with the parameters: 


e v varieties 

e b blocks 

e k varieties in each block, where k < v 
e r occurrences of each variety. 


The variety letters are v and r, and the block letters are b 
and k. 


Theorem 1.1 
In an equi-replicate incomplete block design, the 


parameters v, r, b andk are related by the equation 
or = bk. 


3  Isomorphic designs 

Two block designs are isomorphic if one can be obtained 
from the other by relabelling the blocks and/or the 
varieties. 


4 Connected designs 

Two varieties are directly comparable if there is a block 
that contains both. . 

Two varieties A and B are indirectly comparable if there is 
a chain of varieties A,, ..., A,, such that A and Aj, 
A, and A2,...,and A, and B are all directly comparable. 


A block design is connected if each pair of varieties is 
directly or indirectly comparable. 


5 Incidence matrices 
The incidence matrix B of a block design with v varieties 
and b blocks is the v x b matrix with rows labelled by the 
varieties and columns labelled by the blocks, in which the 
entry in row A and column B is 

1 if variety A lies in block 5; 

Q if variety A does not lie in block B. 
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6 Constructing block designs 

All-combinations method 

For the blocks, take all possible combinations of k varieties 
chosen from the set of v varieties. 


Circle construction 

Place the letters representing the varieties around a circle, 
and mark a starting point; 

for the first block, take the first k varieties, proceeding 
clockwise from the starting point; 

for the second block, take the next k varieties, still 
proceeding clockwise; 

and so on. 


When the varieties are letters, rather than numbers, the 
construction is similar. 


Cyclic construction 

Denote the varieties by 0,1, ..., v-1; 

for the first block, choose any k of these integers; 

for the second block, add 1 to each variety in the first 
block, replacing v by 0 if it occurs; 

for the third block, add 1 to each variety in the second 
block, replacing v by 0 if it occurs; 

and so on. 


When the varieties are letters, rather than numbers, the 
construction is similar. 


7 New designs from old 

Complement of a design a 

From a given design A, the complement A is constructed 
as follows. 3 

For each block 5 of A, construct a block B of A by taking 
the varieties in block B of A to be those varieties that do 
not occur in block B of A. 


The complement of A is A. 


Dual of a design 

From a given design A, the dual A* is constructed as 
follows. 

Interchange the roles of varieties and blocks, re- 
interpreting ‘contains’ as ‘occurs in’, and vice versa. 


The dual of A* is A. 


Section 2 Balanced designs 


1 Concurrence 

Let A and B be two distinct varieties in an incomplete 
block design. The concurrence of A and B is the number of 
blocks where both A and B occur, and is denoted by A,p. 


An incomplete block design A is balanced if there is a 
constant number A such that each pair of distinct varieties 
of A occur together in exactly A blocks. 


The number / is the concurrence of the balanced design A. 


The concurrence matrix C of an incomplete block design is 
the v x v matrix with rows and columns labelled by the 
varieties, in which the entry in row A and column B is: 

r if the varieties A and B are the same; 

Aap if the varieties A and B are different. 


The concurrence matrix of a balanced design with 
replication r and concurrence / has diagonal entries r and 
non-diagonal entries A. 


Theorem 2.1 
The concurrence matrix C of a @regy incomplete block 


design is BB!, where B is the incidence matrix of the 
design. 


If A is any variety, then 
replication of A = number of 1s in row A of B 
= corresponding diagonal entry of BB’ 


If A and B are different varieties, then 
Aap = number of blocks containing both A and B 
= number of 1s in the same position in rows A 
and B of B 
= corresponding non-diagonal entry of matrix BB’ 


2 Balanced incomplete block designs 


Balanced incomplete block design is frequently 
abbreviated to BIBD. 


Theorem 2.2 
Let A be a balanced incomplete block design with v 


varieties, replication r, block size k, and concurrence A. 
Then 


Mov -1) =rk—1). 


The conditions 


k<v, vr=bk, r(k—\)/( —})) isan integer, 


are necessary for a balanced design: they are not sufficient. 


Theorem 2.3: Fisher’s inequality 


In any balanced incomplete block design with v 
varieties and b blocks, 


usp. 


A balanced incomplete block design is symmetric if b = v. 


Theorem 2.4 


For a symmetric BIBD in which v is even, the number 
r — A must be the square of an integer. 


Theorem 2.5 


Let A be a BIBD, and let A be its complement. Then A 
is also a BIBD. 


Theorem 2.6 


Let A be a BIBD, and let A* be its dual. Then A* is a 
BIBD if and only if A is symmetric. 


3 Cyclic balanced designs 


Theorem 2.7 


The rows of the concurrence matrix of a cyclic design 
are obtained by taking the first row and shifting it one, 
two, three, ... places to the right, each time moving the 
right-hand number back to the beginning. 


Arithmetic modulo n is performed on the numbers 
{0, 1, 2, ...,n — 1} by using ordinary arithmetic and then 
adding or subtracting an appropriate multiple of n to 
obtain a number in the set. 


A table of differences is constructed by writing the 
starting block of a design on the left and at the top of the 
table; the entries in the table are obtained by subtracting 
each element at the top from each element on the left, 
using arithmetic modulo v, where v is the number of 
varieties. 


Theorem 2.8 


The first row of the concurrence matrix of a cyclic 
design is obtained by counting the number of times 
each variety occurs in the main body of the table of 
differences. 


Theorem 2.9 


A cyclic design is balanced if and only if the first block 
has the property that all non-zero varieties occur 
equally often among its differences. 


A first block which gives rise to equal frequencies for the 
differences is called a perfect difference set. 


A perfect difference set (modulo v) is a set of distinct 
numbers 

{b, ,b>,...,b4} (modulo v) 
such that the non-zero differences 

b, —b,, b, —b3,.. 
include each non-zero number (modulo v) equally often. 
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Theorem 2.10 
If v is a prime number of the form 4n + 3, where n is a 
non-negative integer, then the set of non-zero squares 
(modulo v) is a perfect difference set. 


4 Constructions for BIBDs 

Construction 1 

Let p be a prime number of the form 4n + 3. 

To obtain a symmetric balanced design with parameters 
v=b=p,- P=k=(p-1)/2,- A= -3)/4 

take as first block the set of non-zero squares (modulo p); 

then construct the other blocks using the cyclic 

construction. 


Construction 2 
Let p be a prime number of the form 4n + 1. 
To obtain a balanced design with parameters 

v=p, b=2p, r=p-1, k=(p-D/2,. A=(-3)/2, 
take the set of non-zero squares (modulo p), and the set of 
non-squares (modulo p), and use the cyclic construction in 
each case; then put the two designs together. 


Construction 3 
Let p be a prime number of the form 4n + 1. 
To obtain a balanced design with parameters 

Qepg+r ld Peze r=. deta) / 2-0 Sipa WZ, 
take the set of squares (modulo p) including 0, and the set 
of non-zero squares (modulo p) together with an extra 
variety z, and use the cyclic construction in each case; then 
put the two designs together. 


Construction 4 
Let p be a prime number of the form 4n + 3. 
To obtain a balanced design with parameters 

papel, b= 2p, tr =p, k=(—p.+)/2,-A=Qp—1/2, 
take the set of non-zero squares (modulo p) together with 
an extra variety z, and the set of non-squares (modulo p) 
together with 0, and use the cyclic construction in each 
case; then put the two designs together. 


Section 3 . Special types of balanced 
design 


1 Steiner triple systems 

A Steiner triple system is a BIBD with k=3 and A= 1. 
In any Steiner triple system: 

°e r=(v-1)/2; 

e v is odd; 

¢ b=v(v -1)/6. 


Theorem 3.1 


In any Steiner triple system, v has the form 6n + 1 or 
6n + 3 for some integer n. 
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2 Finite projective planes 

A finite projective plane is a regular incidence structure 

in which: 

(a) for each pair of distinct points, there is exactly one 
line incident with both; 

(b) for each pair of distinct lines, there is exactly one 
point incident with both. 


A finite projective plane of order n is a symmetric BIBD 
with the parameters: 
v=b=nt+n+1, 


r=k=n+1, A=1. 


Theorem 3.2 


The above two definitions of a projective plane are 
equivalent. 


Theorem 3.3 


If n is a power of a prime number, then there exists a 
finite projective plane of order n. 


3 Latin square designs 

A latin square of side n is an n x n square with n symbols 
arranged in such a way that each symbol occurs just once 
in each row and just once in each column. 


4 Constructions for latin squares 

Cyclic construction 

In the first row, write down n letters in any order. 

In each subsequent row, shift the letters one place to the 
right, moving the last one to the beginning. 


Extended cyclic construction 

In the first row, write down n letters in any order. 

In each subsequent row, shift the letters r places to the 
right, moving the last r letters to the beginning; r can be 
any number from 1 to n such that r and n have no common 
factor. 


Steiner triple system construction 

Take a Steiner triple system with n varieties labelled 
0,1, ...,n- 1; label the rows and columns of the latin 
square similarly, and write 0, 1, ..., 1 -1 down the main 
diagonal. 

For each pair of varieties A and B, find the block of the 
Steiner triple system containing the varieties A and B, and 
place the third variety of that block in row A and column B 
of the latin square. 


5 Uses of latin squares of side n 

1 Two systems each containing n complete blocks, each 
block of one system having just one plot in common 
with each block of the other system, and n varieties. 

2 One system containing n blocks of size n, and n° 
varieties, which are all the possible combinations of 
two treatment factors, each of which has n different 
‘levels’. 


Section 4 Resolvable designs 


1 Resolvability 

A replicate is a set of plots containing each variety exactly 
once. 

An incomplete block design is resolvable if its blocks can 
be grouped into replicates. 


Theorem 4.1 
In a resolvable incomplete block design, in which the 


number of blocks in each replicate is s 
(a) vu=sk; (b) b=rs. 
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Theorem 4.2 


In any resolvable BIBD, b2>v+r-1. 


2 Constructions for resolvable designs 
A (0, 1)-design is a block design in which each con- 
currence is 0 or 1. 


Construction of simple lattice design 

Write down the Kk’ varieties in any order in a k x k square 
array. 

For the first replicate, take the rows of the square as blocks 
of the design. 

For the second replicate, take the columns of the square as 
blocks of the design. 


Construction of triple lattice design 

Construct a simple lattice design with two replicates as 
above. 

Construct a third replicate from the square array as 
follows. 

The first block contains all varieties on the main (top-left 
to bottom-right) diagonal. The other blocks are the other 
‘diagonals’ in the same direction. 


Construction of rectangular lattice design 

Write down the k(k + 1) varieties in any order in a 
(k + 1) x (k + 1) square array with the main diagonal 
missing. 

For the first replicate, take the rows of the square as blocks 
of the design. 

For the second replicate, take the columns of the square as 
blocks of the design. 


3 Orthogonal latin squares and resolvability 

Two n xn latin squares are orthogonal if, when they are 
superimposed, each of the n? possible ordered pairs of 
letters occurs just once. 

Two orthogonal latin squares can be used to obtain 
replicates for the square lattice design. 


To construct a resolvable design with n* varieties 
(numbers) write down an n Xn square array, and two 
orthogonal n x n latin squares (using letters). 

For the first replicate, take the blocks to be the rows of the 
square array. 

For the second replicate, take the blocks to be the columns 
of the square array. 

For the third replicate, associate a letter with each block 
and take the block to be the numbers in the square array 
corresponding to the positions of the chosen letter in the 
first latin square. 

For the fourth replicate, proceed as for the third, using the 
second latin square. 
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Section5 Balanced designs and codes 


Theorem 5.1 


Let A be a BIBD with v varieties, b blocks, replication r, 
block-size k, and concurrence A. Let C be the code 
whose codewords are the rows of the incidence matrix. 
Then the Hamming distance between any two | 
codewords is 2(r — A), and the code C detects up to r—A 
errors and corrects up to r— A-1 errors. 


Index 


[ala + b] construction 34 
across variable 39 
activity network 17 
vertex represents activity 17 
vertex represents event 18 
addition rule 3 
adjacency matrix 
of digraph 7 
of graph 7 
adjacent 5, 6, 7 
algorithm 3 
activity network 
construction 17 
assignment problem 26 
bottleneck assignment 
problem 26 
critical path construction 18 
critical path scheduling 19 
protection scheme 19 
bubble sort 35 
DEPTH-FIRST SEARCH 36 
edge colouring 24 
efficiency of 4 
Eulerian trail in digraph 16 
exhaustion 3 
exponential-time 4 
feasibility 10 
first-fit decreasing 
packing 19 
first-fit packing 19 
Gomory and Hu 15 
greedy 3 
GROW-TREE 36 
heuristic 3 
Hungarian, 
assignment problem 28 
transportation problem 30 
knapsack problem 38 
Kruskal’s greedy 14 
latest starting times 18 
longest path 17 
manipulating queue 37 
manipulating stack 36 
maximum flow 10 
maximum matching 26, 27 
merge two sorted lists 35 
next-fit packing 19 
north neighbour 37 
of Christofides and Viola 4 
polynomial-time 4 
Prim’s greedy 15 
SEARCH 36 
SEARCH 1 35 
SEARCH 2 35 
shortest path 17 
transportation problem 26 
travelling salesman 
problem 38 
vertex colouring 24 
alkane 3 
all-combinations method 44 
AND tree 37 
AND vertex 37 
AND/OR tree 37 
animal (polygonal) 12 
n- 12 
antiprism 13 
arc 3, 7,9 
backward 9 
forward 9 
saturated 9, 10 
unsaturated 9 
arc-disjoint st-paths 8 
Archimedean solid 13 


arcs, multiple 7 

arithmetic modulo n 45 

assignment problem 26, 28 
bottleneck 26 

asymptotically dominate 35 


B 44 
backward arc 9 
balanced binary tree 36 
balanced block design 44 
basic geometric element 11 
basic network 9 
bay 20 
BIBD 45 
bicentral tree 14 
bicentre (of tree) 14 
bin-packing problem 19 
binary digit 3 
binary search tree 36 
binary tree 14 
balanced 36 
height of 36 
binary tree data type 36 
binary trees, counting 14 
binary word 3 
bipartite graph 3, 6 
complete 6 
bit 32 
block 44 
block code 32 
block design 
balanced 44 
connected 44 
cyclic 45 
equi-replicate 44 
incomplete 44 
resolvable 46 
symmetric 45 
block table 13 
Boolean function 34 
Boolean variable 34 
bottleneck assignment 
problem 26 
algorithm 26 
bound 37 
brace 20 
bracing 20 
branch 40 
branch-and-bound method 37 
branch-and-bound tree 38 
branching tree 14 
breadth-first search 37 
spanning tree 37 
breakthrough 10, 27, 28, 30 
bridge 8 
Brooks’ theorem 23 
bubble sort 35 
building-up method for 


constructing spanning tree 14 


C 44 
Cr 40 
capacitance 39 
capacitor 39 
capacity 3 
lower 10 
of arc 9 
of cut 9, 10 
restriction on vertex 9 
upper 10 
capacity rule 15 
Cayley’s theorem 14 
central tree 14 


centre (of tree) 14 
check bit 32 


x’(G) 24 
(Ky) 24 


x%(G) 23 
chord 40 
Christofides and Viola, 
algorithm of 4 
chromatic index 24 
chromatic number 23 
circle construction for block 
design 44 
closed trail 6, 7 
closed walk 6, 7 
co-tree 40 
code 32 
- associated with design 46 
block 32 
cyclic 33 
dual 33 
even-weight 32 
extended 33 
extended Golay 34 
Golay 34 
Hamming 33 
linear 32 
Mariner 9 34 
n-fold repetition 32 
(n, k) 32 
perfect 33, 34 
Reed-Muller 34 
simplex 33 
codeword 32 
column (in framework) 20 
combinatorial problem 3 
combinatorics 3 
compatible edges 23 
complement of design 44 
complement of graph 5 
complete bipartite graph 6 
complete grapIt%, 8, 11 
complete network 15 
component ®& 
multi-terminal 39, 42, 43 
three-terminal 43 
two-terminal 42 
component equation 39 
concurrence matrix 44 


concurrence of block design 44 
concurrence of two varieties 44 


condensing theorem 15 
condensing vertices 15 
conductance 39 
congruent shapes 11 
connected 7 

strongly 7 
connected block design 44 
connectivity 8 

edge 8 

of kinematic pair 21 

optimal 9 

vertex 8 
connector 

maximum 15 

minimum 14 
constraint 21 
construction problem 3 
contraction of graph 22 
convex hull 11 

graph of 11 
convex set 11 
coordinate 11 
cost 3 


47 


cost matrix 26 
counting (trees) 
binary trees 14 
counting problem 3 
country 23 
critical path 18 
construction algorithm 18 
scheduling algorithm 19 
cross polytope 11 
cube 5, 13 
cube graph 11 
cut 9 
capacity of 9 
minimum 9 
cut tree 15 
cut vertex 8 
cutset 8 
cutset rank 40 
cutting-down method for 
constructing spanning tree 14 
cycle 6, 7 
Hamiltonian 6 
cycle graph 5 
cycle method for planarity 
testing 23 
cycle rank 40 
cycle, Hamiltonian 15 
minimum-weight 15 
cyclic block design 45 
cyclic code 33 
cyclic construction for block 
design 44 
cyclic construction for latin 
square 46 
cylindric pair 20 


d 32 
D¢ 40 
d(x, y) 32 
data type 
binary tree 36 
list 35 
stack 35 
decoding 32 
decomposition of graph 25 
into spanning subgraphs 25 
into spanning trees 25 
deg f 22 
degv5 
degree 
of face 13, 22 
of line 13 
of point 13 
of vertex 5, 13 
degree sequence 5 
demand vertex 26 
demi-regular tiling 12 
depth-first search 35, 36 
spanning tree 36 
digraph 3, 7 
adjacency matrix of 7 
connected 7 
disconnected 7 
Eulerian 7 
Hamiltonian 7 
incidence matrix of 7 
simple 7 
strongly connected 7 
sub- 7 
weighted 3 
digraphs 
isomorphic 7 
same 7 
dimension 11 
dimension (of code) 32 


48 


direct graph 20 
directly comparable varieties 44 
disconnected : 
digraph 7 
graph 6 
disjoint sets 3 
divide-and-conquer method 36 
dodecahedron 5, 12 
dom (G) 24 
dominates 35 
dominating set, 
minimum 24, 37 
dual 13 
code 33 
design 44 
incidence structure 13 
network 39 
of connected planar 
graph 22, 23 
polyhedron 12 
tiling 12 
duality theorem 33 
dualize 23 
dummy activity 18 
dummy vertex 26 


earliest starting time 18 
edge 3,5 
edge colouring 24 

greedy algorithm for 24 
edge connectivity M(G) 8 
edge decomposition 25 
edge variable 39 
edge-adjacent polygons 12 
edge-disjoint st-paths 8 
edge-to-edge tiling 11 
edges, multiple 5 
efficiency of algorithm 4 
encoding 32 
enumeration problem 3 
equi-replicate block design 44 
equivalent codes 33 
error syndrome 33 
Euler’s formula for planar graphs 
22 
Euler’s polyhedron formula 13 
Eulerian 

digraph 7, 16 

graph 6, 7 

trail 6, 7, 16 
even-weight code 32 
event 18 
exhaustion algorithm 3 
existence problem 3 
exponential-time algorithm 4 
extended code 33 
extended Golay code 34 


face 22 
degree of 22 
infinite 22 
of polyhedron 12 
face-regular planar graph 22 
factory rule 19 
Fano plane 13 
feasibility 
algorithm 10 
condition 9, 10 
theorem 10 
feasible network 10 
feasible solution 38 
finite projective plane 13, 46 
first order Reed-Muller code 34 
first-fit decreasing packing 
algorithm 19 


first-fit packing algorithm 19 
Fisher’s inequality 45 
five colour theorem 
for maps 23 
for planar graphs 23 
fixed edge 21 
fixed link 20 
float 18 
flow 9, 10 
conservation law 9, 10, 39 
maximum 9, 10 
value of 9, 10 
flow-augmenting path 9, 10 
flow-equivalent network 15 
flow-equivalent tree 15 
formula, chemical 3 
forward arc 9 
four colour theorem 
for maps 23 
for planar graphs 24 
freedoms of link 21 
freedoms of two-link kinematic 
system 21 
fundamental cutset 40 
equation 40 
matrix 40 
fundamental cycle 40 
equation 40 
matrix 40 


Gaussian elimination 41 
generalized max-flow min-cut 
theorem 10 
generator matrix 32 
standard form 32 
generator set 32 
Golay code 34 
Gomory and Hu’s algorithm 15 
graph 3,5 
adjacency matrix of 7 
bipartite 3, 6 
complete 5, 8 
complete bipartite 6 
component of 5 
connected 6, 7 
cycle 5 
disconnected 6 
dual 22 
Eulerian 6 
Hamiltonian 6, 7 
immobile 21 
incidence matrix of 7 
k-colourable 23 | 
k-edge colourable 24 
mobile 21 
non-planar 22 
null 5 
of convex hull 11 
overconstrained 21 
path 6 
Petersen 6 
planar 22 
Platonic 5, 12 
regular 5, 8 
semi-Eulerian 6 
semi-Hamiltonian 6 
simple 5 
sub- 5 
underlying 6 
weighted 3 
graphs 
isomorphic 5 
same 5 
greedy algorithm 3 
for edge colouring 24 
for vertex colouring 24 


Hamiltonian 
cycle 6,7 
digraph 7, 16 
minimum-weight 15 
graph 6, 7 
Hamming code 33 
Hamming distance 32 
Hamming inequality 33 
handshaking dilemma 7 
handshaking lemma 
for graphs 5 
for planar graphs 22 
for polyhedra 13 
heuristic algorithm 3 
for travelling salesman 
problem & 
hexagonal packing 
arrangement 11 
hexagonal tiling 12 
hull, convex 11 
Hungarian algorithm for 
assignment problem 26, 28 
transportation problem 26, 30 


icosahedron 5, 12 
ideal independent source 39 
immobile graph 21 
immobile system 21 
in-degree 7 

sequence 7 
incidence matrix 

of block design 44 

of digraph 7 

of graph 7 

of incidence structure 13 

of oriented graph 41 
incidence relation 13 
incidence structure 13 

dual 13 

regular 13 
incident 

from 7 

to 7 

with 5 
inclusion—exclusion principle 3 
incompatible edges 23 
incomplete block design 44 
ind (G) 24 
indeg v 7 
independence number 24 
independent set 24 

maximum 24 
indirectly comparable 

varieties 44 
inductance 39 
inductor 39 
infeasible network 10 
infinite face 22 
interchange graph 20 
intersection of sets 3 
intractable problem 4 
isomer 3 
isomorphic 

block designs 44 

digraphs 7 

graphs 5 


join 5,7 
joint 20 


k-colourable graph 23 
k-colouring 23 

k-cube 6 

k-edge colourable graph 24 


k-edge colouring 24 

k-screen 37 

kinematic pair 20 

kinematic structure 20 

kinematic system 20 

Kirchhoff’s current law 39 

Kirchhoff’s voltage law 39 

knapsack problem 38 
algorithm 38 

Konig’s theorem 24 

Kruskal’s greedy algorithm 3 
for maximum connector 14 
for minimum connector 14 

Kuratowski’s theorem 22 


latest starting time 18 
latin multiplication 16 
latin square 46 
lattice, point 11 
Laves tiling 12 
left subtree 36 
length (of code) 32 
length of walk 6, 7 
line 
degree of 13 
segment 11 
linear code 32 
linearly dependent equations 40 
linearly independent 
cutsets 40 
cycles 40 
equations 40 
link 20 
list 35 
list data type 35 
listing problem 3 
location problem 4 
longest path 18 
algorithm 17 
loop 5, 7 
lower bound 38 
lower capacity of arc 10 


magnify image 37 
majority-logic decoding 32 
map 23 
Mariner 9 code 34 
marriage condition 25 
marriage theorem 25 
matching 25 
maximum 25 
mathematical model 4 
mathematical modelling 
process 4 
matrix 
adjacency 7 
incidence 7 
transpose of 13 
max-flow min-cut theorem 9 
maximum connector 15 
Prim’s greedy algorithm 
for 15 
maximum flow 9 
algorithm 10 
maximum independent set 24 
maximum matching 25 
algorithm 27 
maximum spanning tree 15 
Menger’s theorems and 
corollaries 8 
merge sort 35 
minimally braced rectangular 
framework 20 
minimum bracing 20 
minimum completion time 18 


minimum connector 3, 14 
Kruskal’s greedy algorithm 
for 14 
Prim’s greedy algorithm 
for 14 
minimum cut 9 
minimum distance 32 
minimum dominating set 24 
minimum spanning tree 14 
mixed network 9 
mobile graph 21 
mobile system 21 
mobility 20 
mobility criterion 21 
for graphs 21 
model 4 
modelling process 4 
modified marriage problem 25 
modified marriage theorem 25 
molecule 3 
motion 20 
multi-terminal 
component 39, 42, 43 
flow problem 15 
multiple arcs 7 
multiple edges 5 
multiplication rule 3 


n-animal 12 
n-cube 11 
n-fold repetition code 32 
(n, k) code 32 
n-octahedron 11 
n-simplex 11 
n-sphere 11 
network 3 

basic 9, 10 

complete 15 

feasible 10 

infeasible 10 

mixed 9 

undirected 9 

with lower and upper 

capacities 10 

next-fit packing algorithm 19 
non-planar graph 22 
north neighbour 37 

algorithm 37 
NP-complete problem 4 
null graph 5 


O(g(n)) 35 
octahedron 5 
n- 11 
off-line algorithm 19 
Ohm’s law 39 
on-line algorithm 19 
open trail 6 
opposites 37 
optimal connectivity 9 
optimization problem 3 
OR tree 37 
OR vertex 37 
order hierarchy 35 
order O(g(n)) 35 
order of magnitude, 
same 35 
order of rotation 12 
Ore’s theorem 7 
oriented graph 39 
original vertex 17 
orthogonal latin squares 46 
orthotope 11 
standard n- 11 
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out-degree 7 

sequence / 
outdeg v 7 
overall parity check 32 
overconstrained graph 21 
overconstrained system 21 


packing 11 
parameter (of design) 44 
parity check equation 33 
parity check matrix 33 
standard form 33 
partial graph 26 
path 6, 7 
critical 18 
flow-augmenting 9 
longest 17, 18 
semi-Hamiltonian 6 
shortest 17 
st-& 
path graph 6 
perfect code 33, 34 
perfect difference set 45 
Petersen graph 6 
planar graph 22 
face-regular 22 
planar pair 20 
planar system 21 
planarity, testing for 23 
plane drawing 22 
face of 22 
infinite face of 22 
plane segment 11 
Platonic graph 5, 12 
Platonic solid 5, 12 
plot 44 
point 11 
degree of 13 
lattice 11 
polygon 11 
polygonal animal 12 
polygons, edge-adjacent 12 
polyhedron 11 
convex 12 
dual 12 
regular 12 
semi-regular 13 
polynomial-time algorithm 4 
polytope 11 
convex 11 
cross 11 
standard 11 
pose 20 
potential difference 
conservation law 39 
precedence relation 17 
Prim’s greedy algorithm 3 
for maximum connector 15 
for minimum connector 14 
tabular (matrix) form 15 
prism 13 
prismatic pair 20 
project clock 19 
rotection scheme 19 
weer SX Dapper ir 
quad tree 37 
queue 36 


R(m) 34 

R(n) 32 

r-regular graph 5 

rate (of code) 32 

reduce image 37 

reduced incidence matrix 41 
reducible configuration 23 


50 


Reed-Muller code 34 
reference direction 39 
reference vertex 41 
reflect image 37 
region of network 4 
regular graph 5, 13 
regular incidence structure 13 
regular polygon 11 
regular polyhedron 12 
regular tiling 12 
replicate 46 
replication 44 
resistance 39 

resistor 39 

resolvable design 46 
Reuleaux pairs 20 
right subtree 36 
rooted tree 14 

rotate image 37 
rotation of order k 12 
rounding 25 

row (in framework) 20 
rules of counting 3 


S(n) 35 
saturated arc 9 
schedule 19 
screw pair 20 
semi-Eulerian graph 6 
semi-Eulerian trail 6 
semi-Hamiltonian graph 6 
semi-Hamiltonian path 6 
semi-regular polyhedron 13 
semi-regular tiling 12 
separate s from t 8 
shadow vertex 17 
Shannon’s theorem 24 
shortest path 17 
algorithm 17 
simple digraph 7 
simple graph 5 
simple lattice design 46 
simplex code 33 
simplex, n- 11 
sink 9 
Stiper- 9 
six colour theorem 
for maps 23 
for planar graphs 23 
solid 5 
Archimedean 13 
Platonic 5, 12 
solution by iteration 35 
solution vector 38 
sorting tree 14 
source 9 
super- 9 
space 
dimension of 11 
segment 11 
space complexity function 35 
spanning tree 14, 40 
maximum 15 
minimum 14 
of oriented graph 40 
spatial system 20 
sphere, n- 11 
spherical pair 20 
spherical system 20 
square packing arrangement 11 
st-path 8 
stack 35 
stack data type 35 
standard cross polytope 11 
standard form of generator 
matrix 32 


standard form of parity check 
matrix 33 

state equation 41 

state space tree 37 

state variable 41 

Steiner triple system 45 
construction for latin 

square 46 

string multiplication 16 

strongly connected 7 

subdigraph 7 

subdivision of graph 22 

subgraph 5 

super-sink 9 

super-source 9 

supply vertex 26 

symmetric block design 45 

systematic encoding rule 32 


T(n) 35 
table of differences 45 
Tellegen’s theorem 41 
terminal 39 
tessellation 11 
tetrahedron 5 
theorems 
Cayley’s 14 
condensing 15 
condition for minimally 
braced rigid framework 20 
condition for rigid braced 
framework 20 
condition for strongly 
connected digraph 16 
describing freedoms of 
link 21 
different capacities 
in C(N) 15 
equivalent definitions 
of tree 14 
first mobility criterion 21 
is connected digraph 
Eulerian 16 
mobility criterion for 
graphs 21 
number of walks in 
digraph 16 
Prim's and Kruskal’s 
algorithms 15 
second mobility criterion 21 
smallest capacity in C(N) 15 
splitting Eulerian digraph 
into cycles 16 
thickness 25 
three-terminal component 39, 43 
through variable 39 
tile 11 
tile the plane 11 
tiling 11 
demi-regular 12 
dual 12 
edge-to-edge 11 
Laves 12 
regular 12 
semi-regular 12 
time complexity function 35 
tractable problem 4 
trail 6, 7 
closed 6, 7 
Eulerian 6, 7 
open 6 
semi-Eulerian 6 
transformation to basic 
network 9 
transformer 39 
transhipment problem 26 


transistor 39 
transportation problem 26, 30 
transpose of matrix 13 
travelling salesman problem ®, 15 
travelling salesman problem, 
algorithm 38 
tree 6 
bicentral 14 
binary 14 
branching 14 
central 14 
cut 15 
equivalent definitions of 14 
flow-equivalent 14 
root of 14 
rooted 14 
sorting 14 
spanning 14 
trees, counting 14 
triangle 6 
triangular packing 
arrangement 11 
triple lattice design 46 
two-terminal component 39, 42 


underlying graph 7 
undirected network 9 
unsaturated arc 9 
upper capacity of arc 10 


value of flow 10 
variable 
across 39 
edge 39 
state 41 
through 39 
variety 44 
vertex 3,5, 7 
connectivity K(G) 8 
cutset 8 
degree of 5 
demand 26 
dummy 26 
in-degree of 7 
original 17 
out-degree of 7 
shadow 17 
supply 26 


vertex type 12 
vertex colouring 23, 24 
greedy algorithm for 24 
vertex decomposition 24 
vertex represents 
activity 17 
event 18 
vertex-disjoint st-paths 8 
vertices, adjacent 5, 7 
Vizing’s theorem 24 
extended version 24 


w(x) 32 
walk 6, 7 
closed 6, 7 
length of 6, 7 
weight 3, 32 
weighted digraph 3 
weighted graph 3 


zero word 32 
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